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ANALYSIS AND CONTROL OF ASYMMETRIC VORTEX FLOWS 
AND SUPERSONIC VORTEX BREAKDOWN 

Osama A. Kandil* 

Accomplishments (Dec. 1, 1991— Nov. 30, 1992) 

The accomplishments which have been achieved in the present year covering the period 
from Dec. 1, 1991 until Nov. 30, 1992 are given. These accomplishments include publications, 
national and international presentations, NASA Research Highlights and presentations, and the 
research group supported under this grant. 

I. Conference Papers, Proceedings and Journal Publications: 

1. Kandil, O. A., Kandil, H. A. and Liu, C. H., “Shock/Vortex Interaction and Vortex Breakdown 
Modes,” IUTAM Symposium on Fluid Dynamics of High Angle of Attack, Paper No. T.1.2 , 
University of Tokyo, Tokyo, Japan, September 13-17, 1992 (Invited), (a copy is enclosed ). 

2. Kandil, O. A., Wong, T. C., Sharaf, H. H. and Liu, C. H„ “Recent Advances in Numerical 
Simulation and Control of Asymmetric Flows Around Slender Bodies,” IUTAM Symposium 
on Fluid Dynamics of High Angle of Attack, Paper No. W.1.2 , University of Tokyo, Tokyo, 
Japan, September 13-17, 1992 (Invited), (a copy is enclosed ). 

3. Wong, T. C., Kandil, O. A., and Liu, C. H., “Computation of Wake Vortex Flows and 
Control of Their Effects on Trailing Wings,” ALAA 4429-92-CP, AIAA Atmospheric Flight 
Mechanics Conference , South Carolina, August 1992, Vol. 1, pp. 280—292, (a copy is 
enclosed ). 

4. Kandil, O. A., Sharaf, H. H. and Liu, C. H., “Active Control of Asymmetric Vortical Flows 
Around Cones Using Injection and Heating,” AIAA 92-4426-CP, AIAA Atmospheric Flight 
Mechanics Conference , South Carolina, August 1992, Vol. 1, pp. 244-253, (a copy is 
enclosed ). 

5. Kandil, O. A., Kandil, H. A. and Liu, C. H., “Critical Effect of Downstream Boundary Con- 
ditions on Vortex Breakdown,” AIAA 92-2601 -CP, 10th Applied Aerodynamics Conference , 
Palo Alto, CA, June 22-24, 1992, pp. 12-26, (a copy is enclosed ). 

6. Kandil, O. A. and Liu, C. H„ “Unsteady Vortex Flows and Flow Control Around Slender 
Bodies and Delta Wings,” Workshop on Supermaneuverability, AFOSR, Lehigh University, 
April 9-10, 1992, pp. 383-417. 

7. Kandil, O. A. and Sharaf. H. H., “Recent Advances in Computational Active Control of 
Asymmetric Vortex Flows, “ Fourth International Conference for Fluid Mechanics , Alexan- 
dria, Egypt, April 28-30, 1992, Vol. I, pp. 237-249. 
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8. Wong, T. C., Kandil, O. A. and Liu, C. H., “Three-Dimensional Computational Study of 
Asymmetric Rows Using Navier-Stokes Equations,” Proceedings of Asian Pacific Confer- 
ence on Computational Mechanics , Hong Kong, December 11—13, 1991, pp. 1365—1371, 
(a copy is enclosed) . 

9. Liu, C. H., Kandil, O. A. and Wong, T. C., “Computational Study for Passive Control of 
Supersonic Asymmetric Vortical Rows around Cones,” Impact Journal of Computing in 
Science and Engineering , Academic Press, Inc., Vol. 4, pp. 80-96, March, 1992, (a copy 
is enclosed) . 

10. Kandil, O. A., Wong, T. C. and Liu, C. H., “Numerical Simulation of Unsteady Asymmetric 
Rows around Cones,” Journal of Ruids and Structures, Academic Press, Vol. 6, pp. 
249-265, February 1992, (a copy is enclosed) . 

11. Kandil, O. A., Wong, T. C. and Liu, C. H., “Prediction of Steady and Unsteady Asymmetric 
Vortical Row Around Cones,” AIAA Journal , Vol. 29, No. 12, December 1991, pp. 1269- 
1278, (a copy is enclosed) . 

12. Kandil, O. A., Wong, T. C., Kandil, H. A. and Liu, C. H., “Thin-Layer and Full Navier-Stokes 
Locally-Conical Asymmetric Solutions,” Accepted for Publication to the ASME Journal of 
Ruids Engineering , December 17, 1991, log. No. 3303-RKA, to appear in April 1993. 

13. Kandil, O. A., Kandil, H. A. and Liu, C. H., “Supersonic Quasi- Axisymmetric Vortex Break- 
down,” Accepted for Publication to the ASME Journal of Ruids Engineering , December 17, 
1991, Log. No. 3302-RKA, to appear in April 1993. 

14. Kandil, O. A., Wong, T. C. and Liu, C. H., “Three-Dimensional Navier-Stokes Asymmetric 
Solutions for Cones and Cone-Cylinder Configurations,” Accepted for Publication to the 
AIAA Journal , November 1991, to appear in December 1993. 

15. Liu, C. H„ Wong, T. C. and Kandil, O. A., “Prediction of Asymmetric Vortical Rows 
Around Slender Bodies Using Navier-Stokes Equations,” Japanese Journal of Ruid Dynamics 
Research , to appear January 1993. 

U. National and International Presentations: 

1. “Unsteady High-a Computational Prediction and Control,” Institute for Space and Astronau- 
tical Sciences (ISAS), Tokyo, Japan , September 17, 1992. (Invited) 

2. “Shock/Vortex Interaction and Vortex-Breakdown Modes,” IUTAM Symposium of Ruid 
Dynamics of High Angle of Attack, Tokyo, Japan , September 13-17, 1992. 

3. “Recent Advances in Numerical Simulation and Control of Asymmetric Rows Around 
Slender Bodies,” IUTAM Symposium of Ruid Dynamics of High Angle of Attack, Tokyo, 
Japan , September 13-17, 1992. 







3 


4. “Computational High-a Aerodynamics,” Presentation given to General John Loh (4 star- 
general) of Langley Air Force Base, ODU, August 31, 1992. 

5. “Computation of Vortex Wake Flows and Control of Their Effects on Trailing Wings, AIAA 
Atmospheric Flight Mechanics Conference, Hilton Head, SC , August 10-12, 1992. 

6. “Active Control of Asymmetric Vortical Flows Around Cones Using Injection and Heating,” 
AIAA Atmospheric Flight Mechanics Conference, Hilton Head, SC , August 10-12, 1992. 

7. “Computation and Control of Vortex Wake Flows,” Directorate Review of Theoretical Flow 
Physics Branch, NASA Langley Research Center , Hampton, VA, July 16, 1992. 

8. “Critical Effects of Downstream Boundary Conditions on Vortex Breakdown,” AIAA 10th 
Applied Aerodynamics Conference, Palo Alto, CA , June 22—25, 1992. 

9. “Vortex- Wake Flows,” Division Review of Theoretical Flow Physics Branch, NASA-Langley 
Research Center , Hampton, VA, June 5, 1992. 

10. “Recent Advances in Computational Active Control of Asymmetric Vortex Flows,” Fourth 
International Conference for Fluid Mechanics , Alexandria, Egypt, April 28-30, 1992. 

11. “Prediction and Control of Unsteady Flows for Supermanuverability,” Taiwan Aeronautical 
Research Center — ODU Symposium , MEM Dept., ODU, April 20-21, 1992. 

12. “Physical Issues and Numerical Simulation of Supersonic Vortex Breakdown,” MEM Dept. 
Televised Seminars, ODU , Norfolk, VA, April 171, 1992. 

13. “Unsteady Vortex Rows and Flow Control Around Slender Bodies and Delta Wings,” 
AFOSR Workshop on Supermanueverability, Lehigh Univ., PA , April 9-10, 1992. Also 
Unsteady Aerodynamics Branch Briefing to the Canadian Air Force , NASA Langley, May 
26, 1992. 

14. “Vortex Research Work” Briefing to AFOSR Program Manager, NASA Langley Research 
Center , Hampton, VA, Feb. 12, 1992. 


HI. Papers Submitted or Accepted for Presentation: 

1. Kandil, O. A., Kandil. H. A. and Liu, C. H., “Three-Dimensional Supersonic Vortex 
Breakdown,” AIAA 93-0526, AIAA 31st Aerospace Sciences Meeting , Reno, Nevada, 
January 11-14, 1993. 

2. Kandil, O. A., Kandil, H. A. and Liu, C. H., “Vortex/Shock Interaction on a 65° -Delta Wing 
in Transonic Flow,” Submitted to the AIAA Fluids and Plasma Dynamics , Orlando, Florida, 
July 6-9, 1993. 
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3. Kandil, O. A., Sharaf, H. H. and Liu, C. H., “Active Control of Asymmetric Conical Flow 
Using Spinning and Oscillation,” submitted to the AIAA Fluids and Plasma Dynamics , 
Orlando, Florida, July 6-9, 1993. 

IV. NASA Research Highlights: 

1. “Three-Dimensional Shock/Vortex Interaction and Vortex Breakdown Modes,” Kandil, O. 
A., Kandil, H. A. and Liu, C. H„ NASA RTOP, July 1992. 

2. “Active and Passive Control of Asymmetric Vortical Flows Around Conical Forebodies, 
Kandil, O. A., Sharaf, H. H., and Liu, C. H., NASA RTOP, July 1992. 

V. Research Group: 

This Principal Investigator is assisted by the following persons in the research group: 

1. Dr. T. C. Wong, Research Associate, MEM Dept., Old Dominion University; Prediction 
and Control of Wake- Vortex Flows. 

2. Mr. Hamdy A. Kandil, Ph.D. Candidate, MEM Dept., Old Dominion University; Internal 
and External Shock/Vortex Interaction and Vortex-Breakdown Modes. 

3. Mr. Hazem H. Sharaf El-Din, Ph.D. Student, MEM Dept., Old Dominion University; Passive 
and Active Control Methods for Asymmetric Flows around Conical Forebodies. 
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Summary 


The problems of asymmetric flow around slender bodies and its control are formulated using 
the unsteady, compressible, thin-layer or full Navier-Stokes equations which are solved using 
an implicit, flux-difference splitting, finite-volume scheme. The problem is numerically 
simulated for both locally-conical and three-dimensional flows. The numerical applications 
include studies of the effects of relative incidence, Mach number and Reynolds number on the 
flow asymmetry. For the control of flow asymmetry, the numerical simulation cover passive 
and active control methods. For the passive control, the effectiveness of vertical fins placed 
in the leeward plane of geometric symmetry and side strakes with different orientations is 
studied. For the active control, the effectiveness of normal and tangential flow injection and 
surface heating and a combination of these methods is studied. 


Introduction 


Flow asymmetry around pointed slender cones develops at critical values of relative incidence 
(ratio of angle of attack to nose semiapex angle) due to short-duration transient disturbances 
or forced disturbances. The origin of the transient disturbances may be a transient side 
slip, an acoustic disturbance, or similar disturbances of short duration. The origin of forced 
disturbances is geometric imperfections in the nose or similar disturbances of permanent 
nature. Flow asymmetry produces side forces, asymmetric lifting forces and corresponding 
yawing, rolling and pitching moments that might be larger than those available by the control 
system of the vehicle. Currently, research efforts are devoted for eliminating or alleviating 
flow asymmetry and its corresponding asymmetric loads. Various methods of passive and 
active control are being studied to learn about their control effectiveness. 

In several recent papers by Kandil, et al. [l]-[4], the unsteady, thin-layer, compressible 
Navier-Stokes equations have been used to simulate steady and unsteady, asymmetric vortex 
flows, including their passive control, around cones with different cross-sectional shapes. 
The emphasis of these papers was extensive computational studies of the parameters which 
influence the asymmetric flow phenomenon and its passive control. Since the computational 
cost associated with the solution of three-dimensional-flow problems at reasonable flow 
resolution is very expensive, all the computational solutions were obtained using a locally- 
conical flow assumption. Such an assumption reduces the problem solution to that on 



two conical planes, which arc in close proximity of each other, and hence it reduces the 
computational cost by an order of magnitude. Moreover, such solutions still provide extensive 
understanding of the flow physics since one can use very fine grids for reasonable flow 
resolution. 

In a later paper, by Kandil, et al. [5], the full Navier-Stokes solutions were compared with 
the thin-layer Navier-Stokes solutions. It was shown that the full Navier-Stokes solutions 
produced thicker ffee-shear layers and more vortex-core resolution as compared with those 
of the thin-layer Navier-Stokes equations. In reference [5], a few tentative three-dimensional 
flow solutions were also presented. 

Substantial research efforts have recently been devoted for eliminating or alleviating flow 
asymmetry and its corresponding side force. In the experimental area, several passive-control 
methods [6]-[8] and active-control methods [9H13] have been investigated. Computational 
simulations have also been used to investigate the effectiveness of several passive-control 
methods [l]-[5] and active-control methods [12], [14], [15]. Various methods of passive 
control were demonstrated in the above references which include the use of vertical fins 
along the leeward plane of geometric symmetry, thin and thick side strakes with different 
orientations, and rotatable forebody tips which have variable cross section (from a circular 
shape at its base to an elliptic shape at its tip). It was shown by Kandil, et al. [4] that 
side-strakes control is more practical than the vertical-fin control since the former was more 
effective over a wide range of angle of attack than the former. Moreover, side-strake control 
provided an additional lifting force. However, the effectiveness of the side-strake control 
terminates at very high angles of attack for the considered strake geometry and flow conditions. 

Various active-control methods have been used which include forebody blowing and movable 
forebody strakes. The forebody blowing methods include forward blowing, normal blowing, 
aft blowing and tangential blowing. The main concept of forebody blowing is to control flow 

separation on the forebody and to create yawing forces and moments which can be utilized 
in controlling the body. 

In this paper, we present samples of simulating asymmetric locally-conical and three- 
dimensional flows around cones. Next, we present samples of simulation for passive control 
using a vertical fin and a side strake. Samples of simulating active control using normal and 
tangential flow injection, surface heating and hybrid methods are also presented. 


Highlights for Formulation and Computational Schemes 


Formulation: The asymmetric-flow problems including their passive and active controls 

are formulated using the conservative form of the unsteady, compressible, thin or full Navier- 


Stokes equations in terms of time-independent, body-conforming coordinates. The equations 
are given in Ref. [5] and hence, they are not repeated here. 

The boundary and initial conditions vary according to the problem under consideration. The 
boundary conditions are explicitly satisfied. In general, they include inflow-outflow conditions 
and solid-boundary conditions. For problems of flow asymmetry, where the flow is solved 
throughout the whole computational domain, periodic boundary conditions are used at the 
plane of geometric symmetry of the problem. 

For the asymmetric flow problems around slender bodies and for supersonic inflow-outflow 
boundary, the Riemann-invariant boundary conditions are used. They require that the inflow 
variables be at the freestream conditions, and the conical shock enclosing the body be captured 
as part of the solution. For supersonic outflow boundary, the Riemann-invariant boundary 
conditions require that all flow variables be extrapolated from the interior cells. On the solid 
boundary, without injection or heating, the no-slip and no-penetration conditions are enforced. 
Moreover, the zero normal-pressure gradient and adiabatic boundary conditions are enforced. 
For the active control problems, the mass-flow rate is specified at the body surface for the 
normal injection control and the temperature distribution is specified at the surface for the 
heating control. For the tangential flow injection, the mass flow rate and velocity profile are 
specified at the lip exit. 

The initial conditions correspond to the uniform flow conditions with ui = U 2 = U 3 = 0 on 
the solid boundary. These conditions are used to obtain the asymmetric flow solution. Next, 
the flow control conditions are enforced and the previously obtained asymmetric solution is 
used for the initial conditions of the active control problem. 

Computational Scheme: The implicit, upwind, flux-difference splitting, finite-volume 

scheme is used to solve the unsteady, compressible, full Navier-Stokes equations. The scheme 
uses the flux-difference splitting scheme of Roe which is based on the solution of the ap- 
proximate Riemann problem. In the Roe scheme, the inviscid flux difference at the interface 
of computational cells is split into two parts; left and right flux differences. The splitting is 
accomplished according to the signs of the eigenvalues of the Roe averaged-Jacobian matrix 
of the inviscid fluxes at the cell interface. The smooth flux limiter is used to eliminate os- 
cillations at locations of large flow gradients. The viscous-and heat-flux terms are linearized 
and the cross-derivative terms are eliminated in the implicit operator. The viscous terms are 
differenced using a second-order accurate central differencing. The resulting difference equa- 
tion is approximately factored and is solved in three sweeps in the £ 2 , and £ 3 directions. 
The computational scheme is coded in the computer program “FTNS3D”. 

For the locally-conical flow solutions, an axial station of xi = 1.0 is selected and the 
components of the flowfield vector are forced to be equal between this axial station and 



another axial station in close proximity to xi = 1.0. This ensures that the flow variables are 
locally independent of the axial direction at xi = 1.0 (Kandil, et al. [5]). 

Computational Applications and Discussion 


Asymmetric Steady Flow (locally-conical solution): Figure 1 shows the residual error 

versus the number of iterations, surface-pressure (SP) coefficient, cross-flow velocity and 
total-pressure-loss (TPL) contour for the solution around a 5°-semiapex circular cone. Two 
computer codes (CFL3D and ICF3D) which solve the thin-layer Navier-Stokes equations arc 
used to validate the asymmetric flow solution. The logarithmic-residual-error curve shows the 
stages through which the solution goes until a stable asymmetric steady solution is obtained. 
A grid of 161x81x2 points in the wrap-around, normal and axial directions, respectively, 
with minimum spacing of 10" 4 , has been used. The computational domain extends around 
the body to 21 r where r is the local radius of the cone. For these critical conditions, the 
asymmetry is developed due to random disturbances; such as the machine round-off-error 
for the CFL3D solution. 

Asymmetric Unsteady Flow (locally-conical solution): Figure 2 shows the results for the 
solution of the flow around the same cone using the same grid, where a = 30°. Here, 
the solution is validated by using the thin-layer and full, Navier-Stokes equations using the 
flux-difference splitting (FDS) scheme and the flux-vector splitting (FVS) scheme. All the 
solutions show asymmetric, periodic flow with vortex shedding. The figures show snapshots 
of TPL contours over a half cycle of periodic response. The periodicity is substantiated by 
the first and last snapshot in each row. It is clearly observed that they are mirror image of 
each other over the half cycle. 

Asymmetric Steady Flow (three-dimensional solutions): Figures 3 and 4 show the asym- 
metric flow results for the flows around a 5° -semiapex cone and a 5°-semiapex cone with 
cylindrical after-body configuration. A grid of 161x81x65 points in the wrap-around, nor- 
mal and axial directions, respectively, with minimum spacing of 10' 6 has been used. For the 
cone solution, the spatial flow asymmetry is qualitatively similar to that of the temporal flow 
asymmetry of the locally-conical flow solution of Fig. 2. For the cone-cylinder configuration 
at the same angle of attack and Mach number as those of the cone case, it is concluded that 
the cylindrical afterbody enhances the flow asymmetry at lower Reynolds number. These 
solutions are obtained using the thin-layer Navier-Stokes equations. 


Passive Control Using Fins and Strakes (locally-conical solutions): Figures 5 and 6 

show the solutions using passive flow controls through a vertical fin placed in the leeward 
plane of geometric symmetry (Fig. 5) or side strakes (Fig. 6). It is concluded that the fin 
height, h, must be at least equal to or greater than the height of the free-shear layers in order 


to yield a symmetric flow. The side-strake control is more practical than the fin control since 
it is more effective for high angles of attack than the fin, and moreover it provides additional 
lifting force. These solutions are obtained using the thin-layer Navier-Stokes equations. 


Active-Control Using Normal Flow Injection: Figure 7 shows the history of the locally- 

conicaly full Navier-Stokes solutions for active control around a 5°-semiapex cone. The 
control is achieved by injecting flow from circumferential ports in the circumferential angle 
range of 6 = ± 67.5°. A variable mass-flow-rate injection of maximum rate of 0.03 is used. 
The mass flow rate is proportional to the difference in the surface pressure between the left 
and right sides of the cone. Figure 8 shows the effectiveness of this flow injection as the 
angle of attack is increased up to 30°, where flow asymmetry develops again. The solution 
is obtained using the full Navier-Stokes solver, FTNS-3D code, on a grid of 241x81x2. 


Hybrid Active Control Using Heating and Injection: Since normal flow injection failed 

to yield asymmetric solutions at angles of attack as high as 30°, hybrid methods of active 
control are investigated. Figure 9 shows the effectiveness of hybrid surface heating and 
variable mass flow normal injection for the same cone as the angle of attack is increased. It 
is seen that this method is very promising at high angles of attack. Here, symmetric surface 
pressure and not symmetric flow is obtained at a = 38° and 42°. The surface temperature is 
taken as T s = 5 and the maximum mass flow rate is 0.05. The solution is obtained using 
the full Navier-Stokes solver, FTNS-3D code, on the same grid as that of Fig. 8. 


Active Control Using Tangential Flow Injection: Figure 10 shows the results for active 

control using injection of flow from side lips. The radius of the lower portion is 1.05 that 
of the upper portion. The maximum mass flow rate is 0.2 and a parabolic velocity profile 
is assumed at the lip exit. The solution is obtained by using the full Navier-Stokes solver, 
FTNS-3D, code and a multi-block scheme to grid the lip-flow exit and the upper and lower 
flow regions. It is seen that this method is effective up to 30° angle of attack. 


Concluding Remarks 


The unsteady, compressible, thin-layer and full Navier-Stokes equations have been used to 
solve for asymmetric steady and unsteady, locally-conical and three-dimensional flows around 
circular cones at high angles of attack. Passive and active control methods have been applied 
to study their effectiveness to yield either a symmetric flow or a symmetric surface pressure 
distribution and hence removing the side forces. Work is underway to use other active control 
methods as well as hybrid passive-active control methods. 
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Thin-layer eqs. (FVS) 



Full Navier-Stokes eqs. (FDS) 

Figure 2. Comparison of unsteady asymmetric now solutions with vortex shedding around a circular cone, 
periodic flow response; a = 30°, Woo = 1-8, H r = ^ = *• 
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Figure 5. Passive control of asymmetric (low around a circular cone using a vertical fin. a = 20°, Mrc = 
1.8, R e = 10 5 , h = 0.5 r and r (r = cone local radius). 



Figure 6. 
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Summary 

Computational simulation and study of shock/vortex interaction and vortex-breakdown modes 
are considered for bound (internal) and unbound (external) flow domains. The problem is 
formulated using the unsteady, compressible, full Navier-Stokes (NS) equations which are 
solved using an implicit, flux-difference splitting, finite-volume scheme. For the bound flow 
domain, a supersonic swirling flow is considered in a configured circular duct and the problem 
is solved for quasi-axisymmetric and three-dimensional flows. For the unbound domain, 
a supersonic swirling flow issued from a nozzle into a uniform supersonic flow of lower 
Mach number is considered for quasi-axisymmetric and three-dimensional flows. The results 
show several modes of breakdown; e.g., no-breakdown, transient single-bubble breakdown, 
transient multi-bubble breakdown, periodic multi-bubble multi-frequency breakdown and 
helical breakdown. 

Introduction 

Longitudinal vortex/transverse shock-wave interactions are typical applications which appear 
in transonic and supersonic flows over a strake-wing configuration at moderate-high angles 
of attack, at a supersonic inlet injesting a vortex and inside a supersonic combustor where 
fuel is injected in a swirling jet to enhance fuel-air mixing [l]-[3]. For the strake-wing 
configuration, vortex breakdown is undesirable since it results in the stall phenomenon, and 
hence its occurrence need to be delayed. On the other hand, vortex breakdown for the other 
two applications is desirable since it enhances mixing and stability of the flame [4]-[5], and 
hence its occurrence need to be controlled for the optimum performance. Unfortunately, the 
literature lacks this type of analysis with the exception of the preliminary work of Liu, Krause 
and Menne [6], Copening and Anderson [7], Delery, et al. [1], Kandil and Kandil [8] and 
Meadows, Kumer and Hussaini [9]. 

The first time-accurate NS solution for a supersonic vortex breakdown was developed by 
the present authors in Ref. [10]. We considered a supersonic quasi-axisymmetric vortex 
flow in a configured circular duct The time-accurate solution of the unsteady, compressible 
NS equations was obtained using an implicit, upwind, flux-difference splitting finite-volume 
scheme. A shock wave has been generated near the duct inlet and unsteady vortex breakdown 
has been predicted behind the shock. The predicted flow was characterized by the evolution, 
convection merging and shedding of vortex breakdown bubbles. The Euler equations were 


also used to solve the same problem. The Euler solution showed larger size and number of 
vortex-breakdown bubbles in comparison with those of the NS solutions. The time-accurate 
solution was carried out for 3,200 times steps which are equivalent to a dimensionless time 
of 16. Only one value of Reynolds number of 10,000 was considered in Ref. [10]. 

In a later paper [11], we expanded our study of this flow using time-accurate computations 
of the NS equations with a fine grid in the shock-vortex interaction region and for longer 
computational times. Several issues were addressed in that study. First, we showed the 
effect of Reynolds number on the temporal evolution and persistence of vortex- breakdown 
bubbles behind the shock. In that stage of computations, the conditions at the downstream 
exit were obtained by extrapolating the components of the flowfield vector from the interior 
cell centers. Although the flow was supersonic over a large portion of the duct exit, subsonic 
flow existed over a small portion of the exit around the duct centerline. Therefore, selected 
flow cases were computed using a Riemann-invariant-type boundary conditions as well as 
other boundary conditions at subsonic points of the duct exit [12]. 

In the present paper, we consider shock/vortex interaction and the resulting vortex breakdown 
modes for quasi-axisymmetric and three-dimensional flows. This study covers bound and 
unbound flow domams. For the bound domain, supersonic swirling flow is considered in a 
configured duct, and for the unbound domain supersonic swirling flow that is issued from a 
nozzle into a uniform supersonic flow of lower Mach number is considered. 

Highlights of Formulation and Computati onal Scheme 

Formulation: The conservative, unsteady, compressible, full Navier-Stokes equations in 

terms of time-independent, body-conformed coordinates ( l , ( 2 and £ 3 are used to solve the 
problem. The equations are given in Ref. [11] and hence they are not presented here. Along 
with these equations, boundary conditions are specified at the computational-domain inlet, 
side wall and downstream exit. The downstream exit boundary conditions will be presented 

and discussed in the next section of the computational results. The initial conditions are also 
presented in the next section. 

Computational Scheme: The computational scheme used to solve the unsteady, compress- 
ible full NS equations is an implicit, upwind, flux-difference splitting, finite-volume scheme. 

It employs the flux-difference splitting scheme of Roe which is based on the solution of the 
approximate one-dimensional Riemann problem in each of the three directions. In the Roe 
scheme, the inviscid flux difference at the interface of a computational cell is split into left 
and right flux differences. The splitting is accomplished according to the signs of the eigen- 
values of the Roe averaged-Jacobian matrix of the inviscid flux at the cell interface. The 
smooth limiter is used to eliminate oscillations in the shock region. The viscous and heat-flux 
terms are differenced using second-order spatially accurate central differencing. The resulting 


difference equation is approximately factored and is solved in three sweeps in the f 1 , £ 2 and 
£ 3 directions. The scheme is coded in the computer program which is called “FTNS3D”. 

The quasi-axisymmetric solutions are obtained using the three-dimensional code by forcing 
the flowfield vector to be equal on two axial planes, which are in close proximity of each other. 

Computational Results and Discussion 


a. Quasi-axisymmetric Vortex Breakdown Modes in a Configured Duct: Figure 1 shows 
an axial plane of a configured circular duct. The design of the duct configuration is intended to 
ensure that the supersonic inflow becomes supersonic at the exit. However, as the computation 
will show, a small portion of the duct exit becomes subsonic at certain times for the specified 
inflow conditions. A grid of 221x51x2 in the axial, radial and circumferential directions, 
respectively, is used for the quasi-axisymmetric flow. For the three-dimensional flow, a grid 
of 221x51x49 in the axial, radial and circumferential directions, respectively, is used. The 
minimum grid length is 0.002. The upstream Mach number is kept at 1.75. 


The initial profile for the tangential velocity is given by 


w 

HZ. 


r 


- exp (-£). 


(i) 


where Uoo = 1-74, r m = 0.2 and k e = 0.1. The maximum swirl ratio /?, is at r = 
0.224 and its value is kept at 0.32. The radial velocity, v, at the initial station is set equal 
to zero and the radial momentum equation is integrated to obtain the initial pressure profile. 
Finally, the density p is obtained from the definition of the speed of sound for the inlet flow. 
With these compatible set of profiles, the computations are carried out accurately in time with 
At = 0.0025. The wall boundary conditions follow the typical Navier-Stokes solid-boundary 
conditions. 


Reynolds number, Re = 4,000: Figure 2 shows snapshots of the streamlines and Mach 

contours for the flow case of Re = 4,000. For this value of Reynolds number, a single 
breakdown bubble is seen at t = 5 and it is convected downstreams as time passes. This 
breakdown bubble is formed during the downstream motion of the inlet shock, which reaches 
its maximum downstream displacement at t = 5. Later on, the shock moves upstream, as 
it is seen at t = 8, while the breakdown bubble is convected in the downstream direction. 
Thereafter, the shock stays stationary at the inlet This swirling flow case shows a transient 
single-bubble breakdown flow. The conditions at the exit are obtained by extrapolation from 
the interior cell centers. 

Reynolds number, Re = 20,000: Figure 3 shows snapshots of the streamlines and Mach 

contours for this case. These snapshots show a vortex breakdown mechanism of evolution, 


convection, merging and shedding of bubbles while the inlet shock is moving downstreams, 
then upstreams and finally downstreams. The inlet shock becomes stationary and no more 
bubbles are developed. This swirling flow case shows a transient multi-bubble breakdown 
flow. 

Reynolds number, Re = 100,000: Figure 4 shows snapshots of the streamlines and Mach 

contours for this case. The downstream boundary conditions are obtained by extrapolating 
all the flow variables from the cell centers at the exit. The streamlines snapshots show multi- 
bubble vortex breakdown evolution, convection, merging and shedding. The time-accurate 
integration was carried out up to t = 200 and the solution showed periodic multi-frequency 
cycles of vortex-breakdown bubbles [21]. An example of the merging of vortex breakdown 
bubbles of same sign of vorticity is shown at t = 17. An example of convection and shedding 
of vortex breakdown bubbles is shown at t = 25. Comparing the streamlines solutions at t 
— 25 and t — 89, it is seen that the solutions are almost the same which conclusively show 
that the breakdown process is periodic. The Mach-contours show the dynamics of inlet shock 
motion. In the time range of t = 3-8, the inlet shock moves upstream toward the inlet and 
its central portion exists outside the inlet section at t = 8. In the time range of t = 8-25, 
the inlet shock moves downstream with corresponding evolution, convection, merging and 
shedding of breakdown bubbles. In the time range of t = 25-45, the inlet shock maintains its 
motion in the downstream direction at a slower rate than before, while another shock, which 
is downstream of the inlet shock, appears and also moves in the downstream direction. The 
evolution, convection and shedding slowly continues until t = 66. In the time range of t = 66- 
78, the downstream shock disappears and a large vortex-breakdown bubble appears and moves 
upstream. This motion of the bubble is accompanied by upstream motion of the inlet shock 
~ 78). Later the inlet shock again moves in the downstream direction and the process is 
repeated. An animation movie has been produced for the total dimensionless time of t = 200. 

Exit Riemann Invariant Condition, p^ = 2POC,: In this case, the back pressure at the 

subsonic points of the duct exit is specified to be p h = 2?,* and the other four variables 
are extrapolated from the interior cell centers. Figure 5 shows snapshots of the streamlines 
and Mach contours of the solution. Comparing the present solution with the solution of the 
previous case (Fig. 4), it is seen that the two solutions are similar with the exception that the 
present solution lags that of the first case in time. The reason behind this behavior is that the 
back pressure pj is larger than that of the first case. Moreover, the Riemann invariant type 
conditions at subsonic points let the downstream effects propagate upstream as time increases. 
The existence of the large back pressure which is felt upstream supports the inlet shock and 
keeps it in the inlet region. 

Three-Dimensional Vortex-Breakdown in a Configured Duct: Figure 6 shows the 

streamlines and Mach contours for the three-dimensional vortex breakdown modes in a con- 


figured duct with the same inflow conditions as those of Fig. 4. At the downstream exit, flow 
conditions arc extrapolated from the interior cell centers. It is observed from the vertical, 
axial planes that the vortex breakdown is a three-dimensional phenomenon. The breakdown 
changes from a two-bubble type (t = 2-5) to a multi-bubble type (t = 7-34). At certain 
times, t = 1 1 and 19, it is observed that the multi-bubble breakdown is followed by a spiral 
breakdown. Beyond t = 46, the breakdown becomes a single -bubble type. 


c. Quasi-axisymmetric Vortex-Breakdown of a Supersonic Flow from a Nozzle: Here, a 

supersonic swirling jet at M } = 3.0, which is issued from a nozzle into a supersonic uniform 
flow of Moo = 2.0, is considered. A grid of 221x51x2 in the axial, radial and tangential 
directions, respectively is considered. The computational domain in an axial plane has the 
dimensions of 7.0 x 3.5 in the axial and radial directions, respectively, where the nozzle exit 
radius r = 1. The freestream Reynolds number is 296,000. The inflow profiles are taken 
from the experimental data of Ref. [2] and they are used as quasi-axisymmetric profiles [12]. 
Figure 7 shows snapshots of the streamlines and Mach contours at selected time steps. The 
streamline figures show multi-bubble breakdown at the early levels of time. These bubbles 
develop due to the shock system existing at the nozzle exit. The shock system consists of 
two oblique shocks; one is weak and the other is strong. For t > 55, the shock system and 
the breakdown bubble move upstreams and remain there with low frequency oscillation. 


d. Three-Dimensional Vortex-Breakdown of a Supersonic Flow from a Nozzle: The flow 
of the previous case is considered for three-dimensional computations using a square-cross 
flow computational domain and a circular cross-flow computational domain. For the first case, 
one-half the square side is 3.5 units and for the second case, the radius is 3.5 units. Three 
types of grid are used for the first case (Gird types 1, 2 and 3) and one type of grid is used for 
the second case (Grid type 4), see Fig. 8. This study shows the dominant effect of the grid 
shape and distribution on the vortex breakdown solutions. While grid type 1 shows two large 
breakdown bubbles. Grids 2, 3 and 4 show substantially smaller breakdown bubbles. More- 
over, the conical shock system ahead of the breakdown bubbles shows substantial difference 
in size, location and shape. This study gives an alarming conclusion to the researchers who 
use substantially coarse grids with a rectangular cross-flow domain and a rectangular grid, 
and still claim capturing of physical vortex-breakdown flows. In the present study, the vortex 
breakdown bubbles of grid type 4 are of the same size order as those of grid type 3, although 
the shock system is somewhat different in shape. 

Concluding Remarks 

Computational simulation and study of shock/vortex interaction and vortex breakdown have 
been considered for internal and external supersonic swirling flow. The time-accurate 
computation for full Navier-Stokes equations is used to produce all the present cases. 


Numerous modes of vortex breakdown flows have been captured. The results show the effects 
of Reynolds number, downstream exit conditions, and grid shape, fineness and distribution. 
Three-dimensional vortex breakdown computations show somewhat different results from the 
quasi-axisymmetric counter parts since the three-dimensional flow provides a flow reliefing 
effect in the circumferential direction. 
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Figure 1. Configured duct grid 


Figure 2. Streamlines and Mach contours for a swirling flow with a transient single-bubble 
breakdown. Moo = 1.75, 0 = 0.32, Re = 4,000. 






t ■ 15 


t = 17 


t* 30 




t = 30 


II u4 


Figure 3. Streamlines for a swirling flow with transient multi-bubble breakdown. Moo 
0 = 0.32, Re = 20,000 

































Figure 6. Mach contours and streamlines for a swirling three-dimensional flow in a configured 
circular duct, M*, = 1.75, Re = 100,000, /3 = 0.32 
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ABSTRACT 

The near-vortex-wake flow of a large aspect-ratio rect- 
angular wing is accurately computed by using the thin- 
layer and full Navier-Stokes (NS) equations. The chord- 
wise section of the wing is a NACA-0012 airfoil and its 
tip is round. The computations have been carried out on 
a fine C-0 grid using an implicit, upwind, flux -difference 
splitting, finite-volume scheme. The thin-layer NS results 
have been obtained with and without flux limiters and 
the full NS results have been obtained without flux lim- 
iters. Flow transition from laminar to turbulent is mim- 
iced by tuming-on the Baldwin-Lomax algebraic model 
at an experimentally prescribed chord-station Location of 
0.05. Comparison of the computed results with each other 
and with the experimental data shows that the full NS 
results give the best resolution of the near-vortex-wake 
flow. Next, the strength of the wing-tip vortex has been 
reduced substantially without reducing the lift coefficient 
by using flow-injection from a slot along a portion of the 
wing tip. The flow injection is directed in the wing plane 
at 45° with the wing-tip chord. 

INTRODUCTION 

Recently, there has been an alarming rate of increase 
of the volume of air civil transport at many airports. With 
this status of air traffic congestion, the hazardous effect 
of wake vortices emanating from large aircraft on light 
aircraft, which pass through these vortices during land- 
ing and take-off operations, is becoming a vital concern 
for operations safety. These vortices are characterized 
with high intensity and turbulence, and may produce high 
rolling moments on trailing aircraft which could exceed 
the available roll control of the trailing aircraft More- 
over, the trailing aircraft, under the influence of these 
vortex trails, could suffer a loss of altitude or climb rate 
in addition to structural damages. These vortices may 
persist up to several miles and for long periods of time 
before their decay, and therefore they play a major role 
in sequencing landing and take-off operations at busy air- 
ports. For example, a minimum safe-separation distance 
between aircraft is dictated depending on the vortex inten- 
sity, wind shear, atmospheric turbulence and temperature 
gradient; among others. 
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The literature shows a few experimental and compu- 
tational investigations that attempt to model and analyze 
wake vortex interaction, merging, decay and their haz- 
ardous effects on trailing aircraft Hallock and Eberle 1 
presented a review of the state of the ait of aircraft wake 
vortices covering the research efforts in the United States 
until the mid-seventies. Experimental wind tunnel and 
airport measurements of the vortex wakes were conducted 
by Dee and Nicholas 1 Harvey and Perry 3 , El-Ramly 4 , 
Wood and McWilliams 5 , Gardoz 6 , Cliffone and Lonzo 7 , 
Olwi and Ghazi*, Liu* and Liu et al 10 . 

Mathematical models and computational schemes 
were developed using inviscid analysis by Chorin and 
Bernard 11 , Hackee and Evans 12 , Yates 13 , Iversen and 
Bernstein 14 , and Rossow 15 . Although an inviscid model 
cannot describe the wake aging including its diffusion, it 
is still capable of producing the wake shape and its dy- 
namics. The mathematical modeling used in the above 
references were based on the use of the point vortex 
method to compute the motion of a finite number of point 
vortices which model the vorticity behind a wing. The 
first three-dimensional inviscid model was introduced by 
Kandil, et. al 16 , where the nonlinear vortex-lattice, which 
was also developed by Kandil, et al 17 , was used to com- 
pute the interference flow between wings and the vortex 
wake hazardous effects. 

Viscous modeling of trailing vortices was first in- 
troduced by McCormick, et al 1 *. Viscous interactions 
in vortex wakes and the effects of background turbu- 
lence, wind shear and ground on two-dimensional vor- 
tex pairs were presented by Bilanin, et al 1 ** ®. Liu and 
his co-workers 21 * 24 studied the interaction, merging and 
decay of vortices in two-dimensional space and of three- 
dimensional vortex filaments. To estimate the effects of 
density stratification, turbulence and Reynolds number on 
vortex wakes, an approximate model was recently devel- 
oped by Greene 23 . Later on, Greene and his co-workers 26 
presented selected results of aircraft vortices which in- 
clude a juncture vortex, a lifting-wing vortex and a wake 
vortex. 

It is crmrliiAvt from the above brief literature sur- 
vey that the problem contains several vortex flow regions 
along with several critical parameters influencing the vor- 
tex wake flows. The vortex wake flow of interacting-wing 
problems can be divided into three regions. The first re- 
gion includes the leading wing and its near-vortex-wake 
flow. The second region includes the interacting vortices 
in the wake. The third region includes the far-wake flow 
along with the trailing- wing flow. In Ref. 27, Kandil, 
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Wong and Liu presented results of the near- vortex- wake 
flow for a large-aspect-ratio rectangular wing using the 
thin-layer, Reynolds averaged, NS equations on a C-O 
grid. Next, a small aspect-ratio rectangular wing was in- 
troduced in the near-wake of the leading wing and the 
problem was recomputed to demonstrate the hazardous 
effects of the vortex wake flow. Two interference cases 
were considered. The first was called the “along-track 
penetration through vortex center” and the second was 
called the "along-track penetration between vortices”, see 
Fig. 1. It was demonstrated that in the first case, the 
trailing wing experienced large rolling moments and in 
the second case the trailing wing experienced loss of lift 
In the present paper, the near-vortex-wake flow is 
recomputed using a finer grid than the grid used in Ref. 
27. Moreover, the thin-layer and full NS equations are 
used along with an upwind scheme, wherein the flux 
limiter of the solver is turned on and off. The goal of this 
part of the paper is to accurately compute the vortex wake 
flow and to study the effect of the numerical dissipation 
of the flux limiter on the computed results. In the second 
part of the paper, alleviation of the strength of the tip 
vortex without reducing the lift force is demonstrated by 
using directed flow injection from a slot along a certain 
length of the wing tip. 

FORMULATION 


Full and Thin-Layer Navier-Stokes Equations 

The conservative form of the dimensionless, unsteady, 
compressible, full Navier-Stokes equations in terms of 
time-independent, body-conformed coordinates t 2 and 
t 3 is given by 

dQ dE m a(Ev), 

dt d£ m dt‘ - 0' m-l-3,s — 1-3 (i) 

where 

= ^ m (x il X 2 ,x 3 ) (2) 

Q = j = j0> put, /? u 2, P u 3, pe]‘ (3) 

E m = inviscid flux 

= J [pu m , + <9i{ m p, pu 2 U m 

+ ^“p. PU 3 U m + d 3 £ m p, (pe + p)U m ]‘ (4) 

(Ev) g = viscous and heat - conduction flux in 
direction 

= j[0> » dkC^AC^i, 

dkrK-rYB - qjt)]*; k = l — 3, n = l — 3 


= 0tf“u k (6) 

The three momentum elements of Eq. (5) are given by 

; j = l - 3 (7) 

The last element of Eq. (5) is given by 


[(«• V 

1 Q „d(a 2 )l 

+ (73ijp;^VJ :p = 1 - 3 (8) 

The single thin-layer approximations of the full Navier- 
Stokes equations demand that we only keep the deriva- 
tives in the normal direction to the body, £ 2 , in the viscous 
and heat flux terms in Eqs. (1), (7) and (8). Thus, we let 

s = 2 for the term a ^» in Eq. (1) and s = 2 and n = 2 
in Eqs. (7) and (8). These equations reduce to 


AQ , 0E, n 

dt dt* at 2 ~ 

+^[^( u I+ u 2 + u j) 
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where 
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( 10 ) 
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= ditAt 2 , y- = w = a p ^ 2 u p (12) 


In Eqs. (1)-(12), the dimensionless variables are refer- 
enced to their appropriate frees tream values. The dimen- 
sionless density p, cartesian velocity components m, u 2 
and U3, total energy per unit mass, e, dynamic viscosity, p 
and speed of sound, a, are defined as the ratio of the cor- 
responding physical quantities to those of the freestrcam; 
namely, p x , a«, p* a^, p x and a,*; respectively. The 
pressure, p, is non-dimensionalized by p^a^., and is re- 
lated to the total energy far a perfect gas by the equation 
of state 


(5) 


P = (7 ~ l)p(* - ^UjUj^ ; j = 1 - 3 (13) 


where 7 is the ratio of specific heats and its value is 1.4. 
The viscosity, /i, is calculated from the Sutherland law 

/i = T3/2 (fTr)’ c = 04317 (14) 

where T is the temperature which is non-dimensionalized 
by Toe- The Prandtl number, P r is fixed at 0.72. The 
Reynolds number is defined as R, = PoeUocL/Moc and 
the characteristic length, L, is chosen as the root chord 
of the wing. 

In Eqs. (1)-(13), the indicial notation is used for 
convenience. The subscripts k, n, p and j are summation 
indices, the superscript or subscript s is a summation 
index and the superscript or subscript m is a free index. 
The partial derivative ^ is referred to by 5k. 

Turbulent Flow 

For the turbulent flow, the Navier-Stokes equations 
are transformed to the Reynolds-averaged equations by 
replacing the coefficient of molecular viscosity, /1, and 
the coefficient of thermal conductivity k with 

nt = a 1 + — f*(i + pi//0 (15) 



where is the effective viscosity, k« the effective ther- 
mal conductivity, im the turbulent viscosity, P r the lami- 
nar Prandtl number, P rt the turbulent Prandd number and 
C p the specific heat under constant pressure. The turbu- 
lent viscosity im is obtained by using the two-layer al- 
gebraic eddy viscosity model which was first developed 
by Cebeci 28 for the boundary-layer equations and modi- 
fied later by Baldwin and Lomax 29 for the Navier-Stokes 
equations. The details of the turbulent model is given in 
Ref. (30) by Wong, Kandil and Liu. 

Boundary and Initial Conditions 

Boundary conditions are explicitly implemented. 
They include the inflow-outflow conditions and the 
solid-boundary conditions. At the plane of geometric 
symmetry, periodic conditions are used. The inflow- 
outflow boundary conditions are implemented using the 
one-dimensional Riemann-invariant conditions normal to 
these boundaries. On the solid boundaries, the no-slip 
and no-penetration conditions are enforced; Ui = 112 = U3 
= 0, and the normal pressure gradient is set equal to zero. 
For the temperature, the adiabatic boundary condition is 
enforced at the solid boundaries. 

The initial conditions correspond to the freestream 
conditions with ui = U2 = U3 = 0 on the solid boundaries. 


The freestream conditions are given by 
Poc = S-OC — Toe = 1, 

Uj,oc ~ Moc 006 0 , 

U 2 oc = Mqc sin q, ^ 17 ^ 

uaoc - 0, 

-1/ 1 , 

Poc - 1/7, + 2 

where a is the angle of attack. 

COMPUTATIONAL SCHEME 

The implicit, upwind, flux-difference splitting, finite- 
volume scheme is used to solve the unsteady, com- 
pressible, thin-layer and Full Reynolds-averaged, Navier- 
Stokes equations. The scheme uses the flux-difference 
splitting scheme of Roe which is based on the solution of 
the approximate Riemann problem. In the Roe scheme, 
the inviscid flux difference at the interface of computa- 
tional cells is split into two parts; left and right flux dif- 
ferences. The splitting is accomplished according to the 
signs of the eigenvalues of the Roe averaged-Jacobian 
matrix of the inviscid fluxes at the cell interface. The 
smooth flux limiter is used to eliminate oscillations at 
locations of large flow gradients. The viscous-and heat- 
flux terms are linearized and the cross-derivative terms 
are eliminated in the implicit operator. The viscous terms 
are differenced using a second-order accurate central dif- 
ferencing. The resulting difference equation is approx- 
imately factored and is solved in three sweeps in the 
and £ 3 directions. The computational scheme is 
coded in the computer program “FTNS-3D" 31 . 

COMPUTATIONAL RESULTS AND DISCUSSIONS 

We consider a rectangular wing of aspect ratio, AR 
= 5.9 and a NACA 0012 chordal section, at an angle of 
attack of o=4.64°. The flow Reynolds number, R. is 
3.2x10* (based on the root-chord length) and the Mach 
number, M« is taken as 0.3 for low speed flows. The 
computational domain consists of a hemispherical bound- 
ary which is followed by a cylindrical boundary. The 
hemisphere center is located at the intersection of the 
wing root-chord and its trailing edge. Its radius is 15 
chords and the cylindrical boundary extends 14 chords in 
the downstream direction beyond the trailing edge. A C- 
O grid of 231x65x65 grid points in the stream wise wrap- 
around, spanwise and normal directions, respectively, is 
generated by using transfinite interpolation. The min- 
imum grid spacings normal to the wing surface at the 
leading-edge and tip regions are chosen to be 10" 5 . The 
grid points close to the trailing edge are clustered up to 5 
chords in the downstream direction. The reason is that for 
the present work ow interest is focused on the resolution 
of the near-vortex-wake flow. Figure 2 shows a typical 
grid for the wing and its near wake. The flow is consid- 
ered to be laminar until the chord station of 0.05 from 
the wing leading edge. Thereafter, the turbulent model is 


turned on to mimic the transition to turbulent flow. This 
location of the transitional flow has been experime ntall y 
determined by Yip and Shubert 32 . 

Thin-Layer and Full Navier-Stokes Solutions 

The flow around the rectangular wing is solved us- 
ing the thin-layer NS equations and the full NS equa- 
tions. The computations are carried out three times using 
the same grid described above and the implicit upwind 
scheme. In the first and second times, the thin-layer NS 
equations are solved using the implicit upwind scheme 
with and without the smooth flux limiter. In the third 
time, the full NS equations are solved using the implicit 
upwind scheme without the smooth flux limiter. Next, 
we compare the results of these three computions with 
each other and with the experimental data of Ref. 32. 

Figure 3 shows a comparison of the history of the 
logarithmic residual error and the lift coefficient of these 
computations. The thin-layer computation with a flux 
limiter (TL-limiter) shows a drop of the residual error of 
about five orders of magnitude after 8000 iteration steps. 
The thin-layer computation without a flux limiter (TL-no 
limiter) shows a drop of the residual error of the same 
order of magnitude as that of the TL-limiter computation 
after 10,800 iteration steps. The full Navier-Stokes com- 
putation without a flux limiter (FNS-no limiter) shows 
a drop of the same order of magnitude as that of the 
TL-limiter computation after 8,900 iteration steps. The 
lift coefficient of the three computations is 0.34090 (TL- 
limiter), 0.35354 (TL-no limiter) and 0.35250 (FNS-no 
limiter). In the next table, we show comparison of the 
three computations and the experimental data for the nor- 
mal force coefficient, Cn, the pitching-moment coefficient 
about the leading edge, C M , and the drag coefficient, C D . 
The results of the case of next section for the tip-jet con- 
trol are also included in the table. The results show that 



Cn 

Cm 

Cd 

TL-limiter 

0.34128 

-0.07581 

0.01857 

TL-no limiter 

0.35374 

-0.08161 

0.01682 

FNS-no limiter 

0.35272 

-0.08137 

0.01702 

Experiment 32 

0.35 

-0.0825 

N/A 

FNS-no limiter 
(tip jet control) 

0.3682 

-0.08667 

0.01732 


the Cn of the FNS-no limiter computation is the closest 
to the experimental C N with an error of +0.78%. The 
corresponding error for the TL-no limiter computation is 
+1.07% and for the TL-limiter computation is -2.49%. 
The error in the computed Cm in comparison with the 
experimental value is —1.37% for the NS-no limiter com- 
putation, -1.08% for the TL-no limiter computation and 
-8.11% for the TL-limiter computation. 

Figure 4 shows a comparison of the computed surface- 
pressure-coefficient distribution in the chordwise direction 
at different spanwise stations with the same experimental 


data. The present computed results are m good agreement 
with those of the experimental data with the exception of 
die peak suction pressure at the wing leading edge. This 
is attributed to the simple modeling of the transitional 
flow at this location, wherein a simple algebraic turbulent 
model is turned on. At the spanwise station of 0.9883 in 
the tip region, it is noticed that the computed Cp usin g 
the TL-no limiter and FNS-no limiter computations are 
in excellent agreement with the experimental data than 
that of the TL-limiter computation, particulary in the 
chordwise direction as of X/C=0.3. 

Figure 5 shows a comparison of the spanwise varia- 
tion of the total pressure-loss contours at chord-stations 
range of X/C=4). 9-5.0 covering the evolution of vortex 
wake. The range of the total pressure-loss contours shows 
that the tip vortex reaches its maximum strength very 
close to the downstream side of the trailing edge. There- 
after, the tip-vortex core expands due to the viscous diffu- 
sion and moves inboards while growing in size. It is also 
observed that the wake thickness increases and moves up- 
wards. The TL-no limiter results of the vortex core and 
wake shear layer show less viscous diffusion and more 
inboard motion than those of the TL-limiter results. The 
FNS-no limiter results show a little better resolution of 
the vortex core and wake shear layer than that of the TL- 
no limiter results. However, the inboard motion of the 
vortex core and wake shear layer is the same as that of 
the TL-no limiter results. Hence, it is concluded that the 
flux limiter is responsible for the vortex-core and wake 
shear-layer diffusions and their small inboard motions. In 
Fig. 6, the particle traces for the wing tip flow and the 
vortex wake flow are shown for the TL-limiter and FNS- 
no limiter computati on s, Again, the TL-limiter results 
show the larger diffusion of the tip vortex than that of the 
FNS-no limiter results. 

Active Injection Control of Tip Vortex 

The hazardous effects of the tip vortex on trailing 
aircraft is very significant as it has been demonstrated in 
Ref. 27. Several ideas exist for passive and active control 
of the tip vortex in order to alleviate its hazardous effects. 
The challenging issue here is to diffuse the tip vortex 
without decreasing the lifting force of the leading aircraft. 
Our first attempt to accomplish these objectives is to use 
active injection from a slot along the wing tip. The slot 
is located on the wing-tip surface and extends from the 
chord station X/O02 to the chord station X/C=0.6. The 
direction of the jet blowing is downstream at 45 s which 
is measured from the x-axis. A constant mass-flow-rate 
coefficient of 0.015 is used for the jet blowing throu gh 
the wing-tip slot. The solution of this case is obtained 
by using the FNS-no limiter computation. Two types 
of initial conditions have been used. In the first, the 
control starts from the uniform flow conditions while in 
the second, the control starts from the FNS-no limiter 
solution of the preceding section. The final steady flow 


solution for both cases showed the same results. Here, 
only the results of the first case are shown. 

Figure 7 shows the history of the residual error and lift 
coefficient for 16,000 iteration steps. The spikes shown 
in the residual error graph are due to the step increases 
of the mass-flow rate from 0.005 to 0.01 and finally to 

0.015. The residual error drops three orders of magnitude 
and the lift coefficient shows a value of 0.3680. The lift 
coefficient is 4.4% higher than that of the FNS-no limiter 
computation without injection. The values of Cn, Cm 
and Cd are included in the table of the preceding section. 
The Cd coefficient shows 1.76% increase than that of the 
FNS-no limiter computation without injection. In Fig. 8, 
the pressure coefficient is shown at the span wise stations 
near the wing tip for comparison with the experimental 
data of the case without injection. It is observed that 
only the tip-section C p near the trailing edge shows a 
slight decrease from that of the experimental data. 

Figure 9 shows the spanwise variation of the total 
pressure-loss contours at chord stations from X/C=0.95 
to X/C=5.0. A comparison of these results with those 
of Fig. 5 reveals the excellent effect of the tip-injection 
control in dissipating the lip vortex. This is also seen in 
Fig. 10 of the particle traces in comparison with those 
of Fig. 6. 

Figure 11 shows a blow-up of the total pressure- 
loss contours at the wing tip covering chord stations 
from X/C=0.2 to X/C=0.9. It is observed that the tip 
vortex existing at X/C=0.2 (at the starting location of the 
slot) is quickly diffused by the blown jet. However, by 
the end of the slot length at X/C=0.6, the tip vortex is 
recovering but with small strength. It is concluded that 
tip-injection control is effective in diffusing the tip vortex 
and moreover it is very important to study the effect of the 
injection-slot length in the future. Other parameters need 
to be included in the future study such as the direction of 
injection, the mass-flow rate and the width of injection. 

CONCLUDING REMARKS 

The thin-layer and full NS equtions have been used 
to accurately solve for the near-vortex-wake flow around 
a laige-aspect-ratio rectangular wing. The flow has been 
computed by using the thin-layer equations with and with- 
out flux limiters and the NS equations without flux lim- 
iters, all on the same grid. The results show the substan- 
tial effects of the flux limiter in diffusing the tip vortex 
and the trailing-edge shear layer and in their small in- 
board motions. The full Navier-Stokes computation with- 
out flux limiters show the best results in comparison with 
the experimental data. Also, it has been demonstrated that 
tip-injection control is substantially effective in diffusing 
the tip vortex without reduction of the wing lift coef- 
ficient Although the drag coefficient slightly increased 
the lift to drag ratio increased. Currently, work is under- 
way to conduct an extensive parametric study on the tip- 
injection control and to study its effect on a trailing wing. 
Passive-control methods are also being investigated. 
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Figure I. Types of interference of Trailing Vortex Wake; (1) cross-track penetration, (2) 
along-track penetration between tip vortices, and (3) along-track penetration 
through tip vortex center. 



Figure 2. Typical computational grid for the wing 
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Figure 3. Iteration history of residual error and lift coefficient. 


F.N.S. -no limiter 





Figure 4. Chord wise surface-pressure distribution at span wise stations for an isolated 
wing and comparison with experimental data. 
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Figure 6. Particle traces for the wing tip and its vortex wake. 
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Figure 7. Iteration history of residual error and lift coefficient for the wing with tip- 
injection control (F.N.S. no limiter). 
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Figure 8. Chord wise surface-pressure distribution at span wise stations for an isolated 
w ing with tip-injection control (F.N.S. no limiter). 
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Figure 11. Span wise variation of total pressure-loss contours at chord wise stations on the 
wing with tip-injection control (F.N.S.-no limiter). 
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ABSTRACT 

The effectiveness of certain active-control methods 
for asymmetric flows around circular cones is investigated 
by using computational solution of the unsteady, com- 
pressible full Navier-Stokes equations. Two main meth- 
ods of active control which include flow injection and 
surface heating are used. For the flow-injection-control 
method, flow injection is used either in the normal di- 
rection to the surface or in the tangential direction to 
the surface. For the surface-heating-control method, the 
temperature of the cone surface is increased. The effec- 
tiveness of a hybrid method of flow control which com- 
bines normal injection with surface heating has also been 
studied. The Navier-Stokes equations, subjected to vari- 
ous surface boundary conditions, are solved by using an 
implicit, upwind, flux -difference splitting, finite-volume 
scheme for locally-conical flow solutions. 

INTRODUCTION 

The problems of prediction, analysis and control of 
asymmetric vortical flows around slender pointed bodies 
are of vital importance to the dynamic stability and con- 
trollability of missiles and fighter aircraft. The onset of 
flow asymmetry occurs when the relative incidence (ratio 
of angle of attack to nose semi-apex angle) of pointed 
forebodies exceeds certain critical values. At these crit- 
ical values of relative incidence, flow asymmetry devel- 
ops due to natural and/or forced dis t ur bances. The origin 
of natural disturbances may be a transient side slip, an 
acoustic disturbance, or similar disturbance of short du- 
ration. The origin of forced dis tu r bances is geometric 
perturbations due to imperfections in the nose geometric 
symmetry or similar disturbances of permanent nature. 
In addition to the relative incidence as one of the influ- 
ential parameters for the onset of flow asymmetry, the 
freestream Mach number, Reynolds number and shape of 
the body-cross sectional area are also important param- 
eters. 

In several recent papers by Kandil et aL M , the un- 
steady, thin-layer, compressible Navier-Stokes equations 
have been used to simulate steady and unsteady, asym- 
metric vortex flows, including their passive control, 
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around cones with different cross-sectional shapes. The 
emphasis of these papers was extensive computational 
studies of the parameters which influence the asymmet- 
ric flow phenomenon and its passive control. Since the 
computational cost associated with the solution of three- 
dimensional-flow problems at reasonable flow resolution 
is very expensive, all the computational solutions were 
obtained u s«"g a locally-conical flow assumption. Such 
an assumpti on reduces the problem solution to that on 
two conical planes, which are in close proximity of each 
other, and hence it reduces the computational cost by an 
order of magnitude. Moreover, such solutions still pro- 
vide extensive understanding of the flow physics since 
one can use very fine grids for reasonable flow resolu- 
tion. These studies showed that asymmetric flow solu- 
tions were unique irrespective of the type of flow distur- 
bance; a random disturbance in the form of a machine 
round-off error or a controlled disturbance in the form 
of a short-duration side-slip disturbance. Unsteady asym- 
metric flow solutions with perfectly periodic vortex shed- 
ding were successfully simulated, and the solutions were 
unique irrespective of the computational scheme used. It 
has also been shown that as the Mach number was in- 
creased, the flow asymmetry was decreased, and as the 
Reynolds number was increased the flow asymmetry was 
increased. Moreover, the cross-sectional shape of the 
cone has been shown to be a very influential parameter 
on the flow asymmetry. Circular sections produced very 
strong flow asymmetry and diamond sections produced 
relatively-weaker flow asymmetry. 

In a later paper, by Kandil et aL s . the full Navier- 
Stokes solu tions were compared with the thin-layer 
Navier-Stokes solutions. It was shown that the full 
Navier-Stokes solutions produced thicker free-shear lay- 
ers and more vortex-core resolution as compared with 
those of the thin-layer Navier-Stokes equations. In Ref. 
5, a few tentative three-dimensional flow solutions were 
also presented. 

Substantial research efforts have recently been de- 
voted for ^liminnring or alleviating flow asymmetry 
and its corresponding side force. In the experimen- 
tal area, several passive-control methods 64 and active- 
control methods*' 13 have been investigated. Computa- 
tional simulations have also been used to investigate the 
effectiveness of several passive-control methods 1 ' and 
active-control methods 1 ^ '*• 13 . Various methods of pas- 
sive control were demonstrated in the above references 
which the use of vertical fins along the leeward 

plane of geometric symmetry, thin and thick side strakes 



with different orientations, and rotatable forebody tips 
which have variable cross section (from a circular shape 
at its base to an elliptic shape at its tip). It was shown by 
Kandil et al. that side-strakes control is more practical 
than the vertical-fin control since the former was more ef- 
fective over a wide range of angle of attack than the latter 
Moreover, side-strake control provided an additional lift- 
ing force. However, the effectiveness of the side-strake 
control terminates at very high angles of attack for the 
considered strake geometry and flow conditions. 

Various active-control methods have been used which 
include forebody blowing and movable forebody strakes 
The forebody blowing methods include forward blowing, 
normal blowing, aft blowing and tangential blowing. The 
main concept of forebody blowing is to control flow 
separation on the forebody and to create yawing forces 
and moments which can be utilized in controlling the 
body. 

In the present paper, we investigate the effectiveness 
of two main methods of active control which include flow 
injection and surface heating. The study of flow-injection 
control covers normal and tangetial injection. Moreover, 
a hybrid method of flow control which combines surface 
heating and normal injection methods is also investigated. 
At this stage of research, the flow control is aiming at 
either rendering the vortical asymmetric flow symmetric 
or rendering the surface-pressure distribution symmetric. 
Computational solution of the unsteady, compressible, 

. Navier-Stokes equations is used for the present work 
with the exception of the tangential injection. For the 
tangential injection, the thin-layer NS equations are used. 
The computational applications are focused on circular- 
section cones, and locally -conical flow assumptions are 
used to substantially reduce the computational cost. 

FORMULATION AND COMPUTATIONAL SCHEME 

Full Navier-Stokes Equations 

The conservative form of the dimensionless, unsteady 
compressible, full Navier-Stokes equations in turns of 
tune-independent, body-conformed coordinates f 1 f 2 and 
<T is given by 


(£*) ■ — v iacous and heat— conduction flux in 
direction 

- qk)]‘; k = 1 - 3, n = 1 - 3 (5) 

U m = 0tf“u k (6) 

The three momentum elements of Eq. (5) are given by 

|r] (7) 

The last element of Eq. (5) is given by 

Wvv, -*)■*£:[(«■*<■ 

■Hie reference parameters for the dimensionless form 

L * and Moc for the 

length, velocity, tune, density and molecular viscosity, 
respectively. The Reynolds number is defined as Re = 
£«VocI<//ioc. where the characteristic length, L, is the 
body length. The pressure, p, is related to the total en- 
ergy per unit mass and density by the gas equation 

P = (7 - l)p Je - ^ ( U J + u 2 + U J)J ( 9 ) 

The viscosity is calculated from the Sutherland law 

04317 (10 > 

and the Prandtl number P r - 0.72. In Eqs. (1)-(10), the 
mdicial notation is used for convenience. 


gQ , &E a a(E T ). 

dt " r d£ m dp — 0;m— 1-3, s = l- 3 (1) 
where 

( m ={ m (x l ,x 2 ,xs) (2) 

- _ q l 

Q - J = jlP>PUt,PU2, puj, />e]‘ (3) 

Em = inviscid flux 

= j[/>IWu,U B + 0i{"p,/m a U B 

+&r , p,pujU B + $rf»p,(p e + p)U B ]» (4) 


Boundary and Initial Conditions 

The boundary and initial conditions vary according 
to the problem under consideration. The boundary con- 
ditions are explicitly satisfied. In general, they include 
inflow-outflow conditions and solid- boundary conditions. 
For problems of flow asymmetry, where the flow is solved 
ttooughout the whole computational domain, periodic 
boundary conditions are used at the plane of geometric 
symmetry of the problem. 

For the asymmetric flow problems around slender 
bodies and for supersonic inflow -outflow boundary, the 
Riemann- in variant boundary conditions are used. They 
require that the inflow variables be at the frees trcam con- 
ditions, and the conical shock enclosing the body be 


captured as part of the solution. For supersonic out- 
flow boundary, the Riemann-invariant boundary condi- 
tions require that all flow variables be extrapolated from 
the interior cells. On the solid boundary, without injec- 
tion or heating, the no-slip and no-penetration conditions 
are enforced. Moreover, the zero normal-pressure gradi- 
ent and adiabatic boundary conditions are enforced. For 
the present active control problems, the mass-flow rate 
is specified at the body surface for the normal injection 
control and the temperature distribution is specified at 
the surface for the heating control. For the tangential 
flow injection the mass flow rate and velocity profile are 
specified at the lip exit (shown in Fig. 10 of next section). 

The initial conditions correspond to the uniform flow 
conditions with ui = U 2 = U 3 = 0 on the solid boundary. 
These conditions are used to obtain the asymmetric flow 
solution. Next, the flow control conditions are enforced 
and the previously obtained asymmetric solution is used 
for the initial conditions of the active control problem. 

Computational Scheme 

The implicit, upwind, flux-difference splitting, finite- 
volume scheme is used to solve the unsteady, compress- 
ible, full Navier-Stokes equations. The scheme uses the 
flux-difference splitting scheme of Roe which is based on 
the solution of the approximate Riemann problem. In the 
Roe scheme, the in viscid flux difference at the interface of 
computational cells is split into two parts; left and right 
flux differences. The splitting is accomplished accord- 
ing to the signs of the eigenvalues of the Roe averaged- 
Jacobian matrix of the in viscid fluxes at the cell interface. 
The smooth flux limiter is used to eliminate oscillations at 
locations of large flow gradients. The viscous-and heat- 
flux terms are linearized and the cross-derivative terms 
are eliminated in the implicit operator. The viscous terms 
are differenced using a second-order accurate central dif- 
ferencing. The resulting difference equation is approxi- 
mately factored and is solved in three sweeps in the ( l , 
£ 2 , and ( 3 directions. The computational scheme is coded 
in the computer program “FTNS3D”. 

For the locally-conical flow solutions, an axial station 
of xi = 1.0 is selected and the components of the flowfield 
vector are forced to be equal between this axial station and 
another axial station in close proximity to x\ = 1.0. This 
ensures that the flow variables are locally independent of 
the axial direction at xi * 1.0 (Kandil et aL s ). 

COMPUTATIONAL APPLICATIONS AND DISCUSSION 

For all the computational applications shown in this 
section, a 5° -semiapex angle circular cone at 1.8 Mach 
number and 10 s Reynolds number is considered. These 
flow conditions were considered earlier (Kandil et al. 4, s ) 
for the same cone for asymmetric-flow prediction and 
passive-control methods. The grid is 241x81x2 points in 
the wrap around, normal directions and axial direction, 
respectively. The grid is generated using a modified 


Joukowski transformation with a minimum grid size of 
lCT 4 in the 4 2 direction at the body surface. For the 
tangential flow injection a multi-block grid has been used 
and it is explained later on. 

Asymmetric Flow, a - 20° 

Figure 1 shows the locally-conical flow solution with- 
out any control around the cone at an angle of attack 
of 20°. The figure shows the total-pressure-loss (TPL) 
contours and the surface-pressure (SP) coefficient. The 
surface pressure is presented versus the angle 0, which 
is measured in the clockwise direction from the leeward 
plane of geometric symmetry. This stable asymmetric 
flow solution is obtained after 6,000 iteration steps. The 
source of flow disturbance at these critical flow condi- 
tions is the truncation error. As it has been shown earlier 
(Kandil et al. 1 ), this solution is unique irrespective of the 
type of source of the flow disturbance. 

Normal Injection Control, 0 = - 215° — ► - 67-5°, 
ill = constant = 0.03, a * 20° 

Next, a constant mass-flow injection of rti = 0.03 is 
applied normal to the cone surface. The circumferential 
range of injection ports extends over 0 - -22.5® — ► -67.5°. 
The solution of the previous asymmetric flow case is 
used as initial conditions along with modified surface 
boundary conditions. Figure 2 shows the TPL contours 
and the SP coefficient from the solution of the controlled 
flow. Although die vortical flow is still asymmetric, as 
it is seen from the TPL contours, the SP curve shows 
a symmetric distribution resulting in a zero side force. 
The TPL contours show that a primary vortex of smaller 
height (in comparison to the case of Fig. 1) is still existing 
on the right side and it is connected through a free- shear 
layer with the left side of the body surface. The injection 
flow on the left side of the body decreases the suction 
pressure on that side making its distribution equal to that 
on the right side. 

Normal Injection Control, 0 = 0° - 67-5°, 

di * constant = 0.03, a * 20° 

In this case, the constant mass-flow injection of ih 
= 0.03 is kept fixed while the circumferential range of 
injection ports is extended to cover the range of 0 « 0® 
— -67.5®. Figure 3 shows the TPL contours and SP 
coefficient of the solution of the controlled flow. The 
vortical flow is still asymmetric but the SP curve shows a 
symmetric distribution, which results in a zero side force. 
It is also noticed that the primary vortex on the right side 
moves to a larger height (in comparison with the cases of 
Figs. 1 and 2) and the left-side free-shear layer also moves 
to the same height level. The SP coefficient curve shows 
less negative pressure coefficient on the leeward side 
in comparison with the case of Fig. 2, which produces 
smaller lifting force. 


Normal Injection Control, 9 » - <7i* — . -f <7i # , 
Variable Mass-Flow Rate, O i 0.03, a * 20° 

In this case, the circumferential range of injection 
ports is extended to cover a symmetric range of 9 * - 
67.5° — » +67.5°. Moreover, the injected mass-flow rate 
is made proportional to the difference in the surface pres- 
sure between corresponding points on the right and left 
side of the body (within 9 ■ - 67.5° -» +67.5°). The maxi- 
mum injected mass-flow rate, itw», which corresponds to 
the maximum pressure difference is restricted to rtu.. = 
0.03. Practically, this control method can be achieved by 
sensing the pressure difference between the right and left 
ports and using it as a feedback control in order to inject 
a mass-flow rate which is proportional to that pressure 
difference. Figure 4 shows the results of this controlled 
flow. The history of the residual error versus the number 
of iterations shows the asymmetric solution up to 6,000 
iterations and the evolution of the symmetric controlled- 
flow solution up to 7,200 iterations. The corresponding 
curve of the history of the side force shows that the side 
force becomes zero at the end of 7,200 iterations. Snap- 
shots of the evolution of the symmetric controlled flow 
solution, in terms of the TPL contours, are given at the 
6,200; 6,400; 6,600; 6,800; 7,000 and 7,200 iteration step. 
The vortical flow is rendered symmetric using the present 
method and the SP coefficient curve shows a symmetric 
distribution. 

Surface Heating Control, 9*9—* -180°, 

T« = 7 Toe, a = 20° 

For the heating control, the left surface temperature 
of the cone, in the circumferential range of 9 • 0° — * 
-180°, has been changed from the adiabatic condition to 
a constant temperature condition of surface temperature, 
T, = 7 Toe. where T* is the freestream temperature. 
Starting from the stable asymmetric flow solution at the 
12,000 iteration step, the left surface temperature is raised 
to T, = 5T« and then it is raised further to T, * 7 T*. 
Figure 5 shows the results of this case which include 
the history of the residual error and the corresponding 
history of the side force along with the TPL contours 
after 24,000 iterations. The side-force curve shows that 
the side force decreases as the temperature increases until 
it vanishes for T,* 7 Toe. Although the side force is zero, 
the vortical flow is still asymmetric. This case shows the 
effect of increasing the gas viscosity with heating in order 
to equalize the surface pressure distribution on the right 
and left sides of the body. This is the first time, we know 
of, that such an active control method has been applied 
to control asymmetric flows. 


Injection Flow Control, Variable Maaa-Flow 
Rate, a - 20° - 30° 

In this case, the effectiveness of the normal injection 
control of the case of Fig. 4 (rtw » 0.03, o ■ 20°, $ * 

- 67.5° — ► 67.3°) is studied for higher angles of attack. 
Figure 6 shows the results of this study which includes 
the history of the residual error versus the number of 
iterations, the corresponding history of the side force and 
the TPL contours along with the SP-coefficient curves for 
a = 22°, 24°, 26°, 28° and 30°. The residual error and the 
side-force curves show the responses for the following 
history: Fust, no control is applied until the 12,000 
iteration step; second, injection is applied at a = 20° 
for 2,000 iterations; third, the angle of attack is increased 
2° each 300 iterations until a * 28°; fourth, at a = 30° 
injection is continued for 2,000 iterations. The side-force 
curve shows zero value in the angle of attack range of o = 
20° - 28° (corresponding to the iterations range of 14,000 

- 16,000). When the angle of attack reaches 30°, the 
normal injection control becomes incapable of achieving 
flow symmetry. It is interesting to study the snapshots of 
the TPL contours at these angles of attack. It is noticed 
that the primary vortices increase in size in the normal 
direction and their inner boundaries approach each other. 
At a a 28°, the inner boundaries of the primary vortices 
touch each other, and thereafter at o « 30° the primary 
vortices become asymmetric. The asymmetric response 
of the primary vortices at o - 30® is believed to be due 
to a strong in stab i lit y arising due to the strong interaction 
of the inner shear layers of the primary vortices. 

At a > 30°, the injected mass-flow rate in the normal 
direction to the body surface is increased to tflm,. = 0.03, 
0.06 and 0.07 in order to recover the flow asymmetry. 
Figure 7 shows the results of these tests. Both the TPL 
contours and the SP-coefficient curves show that the flow 
is asymmetric and the side force does not vanish. It is 
seen that the flow asymmetry changes from the right to the 
left sides and so does the side force. Therefore, it is ap- 
parent that normal-injection control loses its effectiveness 
at high angles of attack for the considered conditions. 

Hybrid Heating- Injection Control, a * 20® - 42° 

In this case, the surface-heating control is followed 
by normal-injection control. Figure 8 shows the results 
of this study. The control process is applied as follows: 
Fust, the whole surface is heated to T, = 5T,* and the 
beating control started at a * 20°, where the number 
of iterations is 12,000. Next, the angle of anrir is 
increased 4* until symmetric flow is recovered at the 
18,000 iterations. The angle of attack is then increased 
2 s each 1000 iteration until the angle of attack reaches 
38°. The flow symmetry has been continuously recovered 
until the angle of attack of 36° where a slight flow 
asymmetry appears in the solution. At a ■ 38 s , the flow 
asymmetry becomes stronger and the heating control loses 
its effectiveness for T, = 5Toc. 


At a = 38°, normal-injection control is applied with 

rti = 0.05 and 9 = - 67.5® — 67.5°. The angle of 

attack is then increased 2° each 3000 iterations. The 
history of the residual error and the corresponding zero 
side force is shown in Fig. 8. Figure 9 shows snapshots 
of the TPL contours and SP-coefficient curves at o = 38°, 
40° and 42°. The TPL contours show asymmetric flows 
but the SP-coefficient curves show symmetry resulting 
in zero side force. Thus, it is concluded that a hybrid 
control of heating and normal injection is very effective 
in eliminating the side force at very high angles of attack. 

Tangential-Injection Control, Abu = 0.2 

In this case, the circular section of the lower portion 
of the cone has a raidus larger than that of the upper 
portion. The lower local radius is 1.05 that of the upper 
one. Flow is injected from the lip of the lower portion in a 
direction tangent to the surface of the upper portion. The 
maximum mass flow rate, itw. is 0.2 and the velocity 
profile at the lip exit is a parabolic one. Figure 10 shows 
the grid used for the thin-layer computations. Three 
grid blocks are used for the present computation. The 
first block is 121x19x2, the second is 121x63x2 and the 
third is 121x63x2 all in the circumferential, normal and 
axial directions, respectively. The tangential flow control 
starts from the 20° angle of attack. The computations are 
carried out until a symmetric flow solution is obtained. 
Thereafter, the angle of attack is increased to 30° with 
a step of 2°. Figure 10 also shows the total pressure- 
loss contours for the condoled, symmetric flow solutions 
at o = 22°, 26° and 30®. It is noticed that as the 
angle of attack increases the inner shear layers of the 
condoled symmetric primary vortices approach each other 
and hence flow asymmetry might develop again. 

CONCLUDING REMARKS 

The computational solution of the unsteady, com- 
pressible, full Navier-Stokes equations has been used to 
study the effectiveness of two active control methods for 
asymmetric flows around circular pointed cones. Locally- 
conical flow assumptions have been used in order to re- 
duce the computational cost of this study by an order of 
magnitude in comparison with that of three-dimensional 
flow solutions. The first active control method is applied 
by injecting air normal to the body surface from ports that 
are circumferentially distributed within a certain angle. 
The injected mass-flow rate is either constant or variable. 
The variable injected mass-flow rate is proportional to the 
surface-pressure difference between corresponding points 
on the right and left sides of the body surface. The results 
of this part of the study show that use of normal injection 
alone loses its effectiveness at very high angles of attack. 
The second active control method is applied by heating 
the body surface either partially or totally. The results of 
this part of the study show that surface-heating control is 
much more effective in comparison with normal-injection 


control at very high angles of attack. Next, a hybrid con- 
trol method which includes surface-heating control fol- 
lowed by normal-injection control is applied. The results 
of this part of the study shows the superior control ef- 
fectiveness of the hybrid method. Active control using 
tangential flow injection is very promising. Work is still 
underway to exploit its effectiveness. Simultaneous ap- 
plication of passive control and the present active control 
methods is also investigated. 
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Figure 1. Asymmetric flow solution, TPL contours and SP 
coefficient; a = 20°, = 1.8, Re-10 5 . 


1/ n = 6,200 


120 180 240 300 360 

0 


Figure 2. Normal-injection control, TPL contours and SP 
coefficient; a = 20°, = 1.8, R« = 10 5 , constant 

rti = 0.03, 6 = - 22.5° - 67.5°. 
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Figure 3. Normal-injection control, TPL contours and SP 
coefficient; a = 20°, M<» = 1.8, R* = 10 5 , const rti 
= 0.03, $ = 0° 
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Figure 4. Normal-injection control; history of residual error 
and side force, TPL contours and SP coefficient; 
a = 20°, Moo = 1.8, Re = 10 5 , variable ih, riw = 
- 67.5° — +67.5°. 
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ABSTRACT 

The unsteady, compressible, full Navier-Stokes (NS) 
equations are used to study the critical effects of the 
downstream boundary conditions on the supersonic vortex 
breakdown. The present study is applied to two super- 
sonic vortex breakdown cases. In the first case, quasi- 
axisymmetric supersonic swirling flow is considered in a 
configured circular duel, and in the second case, quasi- 
ax isymmetric supersonic swirling jet, that is issued from 
a nozzle into a supersonic jet of lower Mach number, is 
considered. For the configured duct flow, four different 
types of downstream boundary conditions are used, and 
for the swirling jet flow from the nozzle, two types of 
downstream boundary conditions are used. The solutions 
are lime accurate which are obtained using an implicit, 
upwind, flux-difference splitting, finite-volume scheme. 

INTRODUCTION 

Numerical simulation of vortex breakdown has been 
focused in most of the existing literature on incompress- 
ible flows. Quasi-axisymmctric, vortex-breakdown sim- 
ulations of incompressible swirling flows have been pre- 
sented by Grabowski and Berger 1 , Hafez, et al. 2 , Salas 
and Kuruvilla 3 , Menne 4 , Wu and Hwang 5 and Menne and 
Liu 6 . Three-dimensional, vortex-breakdown simulations 
of incompressible swirling flows have been presented by 
Spall, et al. 7 , Hanel 8 , and Krause 9,10 . Discussion of the ef- 
fects of side boundary conditions and downstream bound- 
ary conditions on vortex breakdown was presented by 
Krause in Ref. 10. 

Longitudinal vortex/transverse shock-wave interac- 
tions are typical applications which appear in transonic 
and supersonic flows over a strake-wing configuration 
at moderate-high angles of attack, at a supersonic in- 
let injesting a vortex and inside a supersonic combustor 
where fuel is injected in a swirling jet to enhance fuel- 
air mixing 11 ' 13 . For the strake-wing configuration, vortex 
breakdown is undesirable since it results in the stall phe- 
nomenon, and hence its occurrence need to be delayed. 
On the other hand, vortex breakdown for the other two 
applications is desirable since it enhances mixing and sta- 
bility of the flame 14 * 15 , and hence its occurrence need to be 
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controlled for the optimum performance. Unfortunately, 
the literature lacks this type of analysis with the excep- 
tion of the preliminary work of Liu, Krause and Menne 16 , 
Copening and Anderson 17 , Delery, et al, n Kandil and 
Kandil 18 and Meadows, Kumer and Hussaini 19 . 

The first time-accurate NS solution for a supersonic 
vortex breakdown was developed by the present au- 
thors in Ref. 20. We considered a supersonic quasi- 
axisymmetric vortex flow in a configured circular duct. 
The time-accurate solution of the unsteady, compress- 
ible NS equations was obtained using an implicit, up- 
wind, flux-difference splitting finite-volume scheme. A 
shock wave has been generated near the duct inlet and 
unsteady vortex breakdown has been predicted behind 
the shock. The predicted flow was characterized by the 
evolution, convection and shedding of vortex breakdown 
bubbles. The Euler equations were also used to solve the 
same problem. The Euler solution showed larger size and 
number of vortex-breakdown bubbles in comparison with 
those of the NS solutions. The time-accurate solution was 
carried out for 3,200 times steps which are equivalent to 
a dimensionless time of 16. Only one value of Reynolds 
number of 10,000 was considered in Ref. 20. 

In a later paper 21 , we expanded our study of this flow 
using time-accurate computations of the NS equations 
with a fine grid in the shock-vortex interaction region 
and for long computational times. Several issues were 
addressed in that study. First, we showed the effect of 
Reynolds number on the temporal evolution and persis- 
tence of vortex-breakdown bubbles behind the shock. In 
that stage of computations, the conditions at the down- 
stream exit were obtained by extrapolating the compo- 
nents of the flowfield vector from the interior cell centers. 
Although the flow was supersonic over a large portion of 
the duct exit, subsonic flow existed over a small portion 
of the exit around the duct centerline. Therefore, selected 
flow cases were computed using a Riemann-invari ant- 
type boundary conditions at subsonic points of the duct 
exit Finally, the effect of swirl ratio at the duct inlet 
was investigated. 

In the present paper, we address the problem of spec- 
ifying the downstream boundary conditions and their crit- 
ical effects on the supersonic vortex breakdown problem 
for internal and external flows. For this purpose, the un- 
steady, compressible, full NS equations are used along 
with an implicit, upwind, flux-difference splitting, finite- 
volume scheme for the time-accurate solutions. For the 


internal flow case, supersonic swirling flow in a config- 
ured duct is considered along with four types of down- 
stream boundary conditions. Keeping the duct geometry 
and the upstream flow conditions fixed, the exit boundary 
conditions are varied. The four exit boundary conditions 
include extrapolation of all the five variables from the 
interior cell centers, specifying the downstream pressure 
and extrapolating the other flow conditions from the in- 
terior cell centers, specifying the downstream pressure 
gradient and extrapolating the other four conditions from 
the interior cell centers, and using a disk of specified 
radius at the exit section. For the external flow case, 
supersonic swirling jet is issued from a nozzle into a su- 
personic non-swirling jet of Mach number lower than that 
of the swirling jet. Two types of downstream boundary 
conditions are considered. In the first type, extrapolation 
of all the five variables from the interior cell centers is 
used, while in the second type, the standard Riemann- 
invariant type boundary condition is used. All the results 
presented in this paper have been produced using our full 
NS solver which is known as the “FTNS3D” code. The 
CRAY-YMP computer of NASA Langley Research Cen- 
ter is used for the computations. 

HIGHLIGHTS OF THE FORMULATION 
AND COMPUTATIONAL SCHEME 

The conservative, unsteady, compressible, full Navier- 
Stokcs equations in terms of time-independent, body- 
conformed coordinates £ 1 , and £ 3 are used to solve the 

problem. The equations are given in Ref. 21 and hence 
they are not shown here. Along with these equations, 
boundary conditions are specified at the computational- 
domain inlet, side wall and downstream exit. The down- 
stream exit boundary conditions will be presented in the 
next section of the computational results. The initial con- 
ditions will also be presented in the next section. 

The computational scheme used to solve the unsteady, 
compressible full NS equations is an implicit, upwind, 
flux-difference splitting, finite- volume scheme. It em- 
ploys the flux -difference splitting scheme of Roe which is 
based on the solution of the approximate one-dimensional 
Ricmann problem in each of the three directions. In the 
Roe scheme, the inviscid flux difference at the interface 
of a computational cell is split into left and right flux dif- 
ferences. The splitting is accomplished according to the 
signs of the eigenvalues of the Roe averaged-Jacobian 
matrix of the inviscid flux at the cell interface. The 
smooth limiter is used to eliminate oscillations in the 
shock region. The viscous and heat-flux terms are lin- 
earized and the cross-derivative terms of the viscous Ja- 
cobians are dropped in the implicit operator. These terms 
arc differenced using second-order spatially accurate cen- 
tral differencing. The resulting difference equation is ap- 
proximately factored and is solved in three sweeps in the 
£ 2 and £ 3 directions. The scheme is used for third- 
order spatial accuracy and first-order temporal accuracy. 


The scheme is coded in the computer program which is 
called “FTNS3D". 

The quasi-axisymmetric solutions are obtained using 
the three-dimensional code by forcing the flowfield vector 
to be equal on two axial planes, which are in close 
proximity of each other. 


COMPUTATIONAL RESULTS AND DISCUSSION 


I. Vortex Breakdown in Configured Circular Duct 

Figure 1 shows a configured circular duct which con- 
sists of a short, straight cylindrical part at the inlet which 
is followed by a short divergent cylindrical part until the 
axial length of 0.74. The divergence angle is 6°. The duct 
radius is then kept constant and a convergent-divergent 
nozzle with a throat radius of 0.95 is attached. The duct 
exit radius is 0.98 and its total length is 2.9. The diver- 
gent part of the duct ensures the stability of the formed 
shock in the inlet region. The configuration of the duct is 
intended to ensure that the supersonic inflow will become 
supersonic at the exit. As the computations will show, a 
small portion of the duct exit near its centerline becomes 
subsonic at certain times for the specified inflow condi- 
tions. This configured duct has also been used by Delery, 
ct al. 11 for their Euler equations computations of super- 
sonic vortex breakdown in an attempt to computationally 
model an experimental set up. 

The NS solver is used to compute all the following 
flow cases by using a grid of 221 x51 on two axial planes, 
where 221 points are in the axial direction and 51 points 
are in the radial direction. In the inlet region up to 
the 0.74 axial station, 100 grid points are used and the 
remaining 121 points are used in the remaining part of the 
duct. The grid is also clustered at the centerline (CL) and 
the wall. The minimum radial grid size at the CL is 0.002. 
The two axial planes are spaced circumferentially at a 
certain angle so that the aspect ratio of the minimum grid 
size will be less than 2. The upstream Mach number is 
kept at 1.75 and the Reynolds number is kept at 100,000. 
The initial profile for the tangential velocity is given by 
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where f/* = 1.74, r m = 0.2 and k t = 0.1. The max- 
imum swirl ratio /?, is at r * 0.224 and its value is 
kept at (T.32. The radial velocity, v, at the initial station 
is set equal to zero and the radial momentum equation is 
integrated to obtain the initial pressure profile. Finally, 
the density p is obtained from the definition of the speed 
of sound for the inlet flow. With these compatible set 
of profiles, the computations are carried out accurately in 
time with At = 0.0025. The wall boundary conditions fol- 
low the typical Navier-Stokes solid-boundary conditions. 
These computations have been carried out on the CRAY 
YMP of the NASA Langley Research Center. The CPU 
lime is 40 /is/grid point/iteration for the NS calculation. 


Next, we present the results of the computational 
study of this case which covers four types of exit bound- 
ary conditions. 

1.1. Extrapolation from Interior Cel! Centers 

Figure 2 shows snapshots of the streamlines and Mach 
contours for this case. The downstream boundary condi- 
tions are obtained by extrapolating all the flow variables 
from the cell centers at the exit. The streamlines snap- 
shots show multi-bubble vortex breakdown evolution, 
convection, merging and shedding. The time-accurate 
integration was carried out up to t = 200 and the so- 
lution showed periodic multi-frequency cycles of vortex- 
breakdown bubbles 21 . An example of the merging of 
vortex breakdown bubbles of same sign of vorticity is 
shown at t = 17. An example of convection and shed- 
ding of vortex breakdown bubbles is shown at t = 25. 
Comparing the streamlines solutions at t = 25 and t = 
89, it is seen that the solutions are almost the same which 
conclusively show that the breakdown process is periodic. 
The Mach-contours show the dynamics of inlet shock mo- 
tion. In the time range of t = 3-8, the inlet shock moves 
upstream toward the inlet and its central portion exists 
outside the inlet section at t = 8. In the time range of 
t = g-25, the inlet shock moves downstream with corre- 
sponding evolution, convection, merging and shedding of 
breakdown bubbles. In the time range of t = 25-45, the 
inlet shock maintains its motion in the downstream di- 
rection at a slower rate than before, while another shock, 
which is downstream of the inlet shock, appears and also 
moves in the downstream direction. The evolution, con- 
vection and shedding slowly continues until t = 66. In 
the time range of t = 66-78, the downstream shock disap- 
pears and a large vortex-breakdown bubble appears and 
moves upstream. This motion of the bubble is accom- 
panied by upstream motion of the inlet shock (t = 78). 
Later the inlet shock again moves in the downstream di- 
rection and the process is repeated. An animation movie 
has been produced for the total dimensionless time of t 
= 200. Figure 3 shows a blow-up of two snapshots of 
the streamlines solutions at t = 84 and t = 87. At t = 84, 
we recognize five vortex breakdown bubbles which spa- 
tially alternate in their sign of vortex strength. It is also 
noticed that six stagnation points exist at the axis. At t 
= 87, we recognize seven vortex breakdown bubbles and 
seven stagnation points. The figure shows the merging of 
two bubbles of same vorticity sign. 

1.2. Riemann Invariant Boundary Conditions 

with p b = poc 

In this case, the back pressure at the subsonic points 
of the duct exist, p bf is specified to be equal to p<x and 
the other four variables are extrapolated form the inte- 
rior cell centers. The computations have repeated on the 
same grid and for the same flow conditions as that of the 
previous case. Figure 4 shows snapshots of the stream- 
lines and Mach contours of the solution. Comparing the 


present solution with the previous case (Fig. 2), it is seen 
that the two solutions are the same until t = 35. There- 
after, for t > 35, the inlet shock continuously moves in the 
downstream direction with the vortex-breakdown bubbles 
moving ahead of the shock. The shock and vortex bubbles 
are shed and disappears from the duct at advanced levels 
of time. The breakdown mode is termed as "a transient 
multi-bubble vortex breakdown." The reason behind dis- 
appearance of the shock-vortex-breakdown-bubble sys- 
tem is that the back pressure is so low that it cannot sup- 
port the inlet shock keeping it in the inlet region. More- 
over, the Riemann-invariant type conditions at subsonic 
points let the downstream effects propagate upstream as 
time increases. 

1.3. Riemann Invariant Boundary Conditions 

with p = 2poc 

In this case, the back pressure at the subsonic points 
of the duct exit is specified to be p b = 2poc and the 
other four variables are extrapolated from the interior cell 
centers. Figure 5 shows snapshots of the streamlines and 
Mach contours of the solution. Comparing the present 
solution with the solution of the first case (Fig. 2), it is 
seen that the two solutions are similar with the exception 
that the present solution lags that of the first case in 
time. The reason behind this behavior is that the back 
pressure p b is larger than that of the first case. Moreover, 
the Riemann invariant type conditions at subsonic points 
let the downstream effects propagate upstream as time 
increases. The existence of the large back pressure which 
is felt upstream supports the inlet shock and keeps it in 
the inlet region. 

1.4. Riemann Invariant Boundary Conditions 

with || = constant 

In this case, the back pressure is obtained from the 
condition || = constant at the subsonic points of the duct 
exit The other flow variables are extrapolated from the 
interior cell centers. Figure 6 shows snapshots of the 
streamlines and Mach contours of the solution. Compar- 
ing the present solution with the solution of the first case 
(Fig. 2), it is seen that the two solutions are similar until 
t = 22. Thereafter, for t > 22, the inlet shock continu- 
ously moves in the downstream direction with the vortex- 
breakdown bubbles moving ahead of the shock. Again, 
as in the case of Fig. 4, the shock and vortex bubbles 
are shed and disappear from the duct at advanced levels 
of time. The breakdown is termed as “a transient multi- 
bubble vortex breakdown." The reason behind disappear- 
ance of the shock-vortex-breakdown-bubble system is that 
the back pressure obtained from = constant condition 
is so low that it cannot support the inlet shock keeping it 
in the inlet region. Moreover, the Riemann-invariant type 
conditions at subsonic points let the downstream effects 
propagate upstream as time increases. 



I. 5. Placing a Disk at tbe Exit with r = 0.333 

In this case, a circular solid disk of radius r = 0.333 is 
placed at the duct exit and solid-boundary conditions are 
applied on the disk surface. For the remaining portion 
of the exit, the boundary conditions are obtained by 
extrapolation from the interior cell centers. Figure 7 
shows snapshots of the streamlines and Mach contours 
of the solution. It is noticed from the Mach contours 
that most of the exit points are subsonic. The streamlines 
show that the shock-vortex-breakown-bubble system first 
appears behind the duct exit. Thereafter, they move 
upstreams until the inlet shock moves outside of the inlet 
The vortex-breakdown bubbles are then locked between 
the shock outside of the inlet and the circular disk at the 
exit 

II. Supersonic Swirling Jet from a Nozzle 

Figure 8 shows the computational domain and a typ- 
ical grid for this external flow case. The dimensions of 
the computational domain is 7 x3.5 in the axial and radial 
directions, respectively. The grid is 221x51x2 points in 
the axial, radial and tangential directions, respectively. 
The grid is clustered at the nozzle exit (x = 0.0) and at 
the domain centerline. The dimensionless nozzle radius 
is unity, where a supersonic jet of M x = 3.0 is issued, 
and outside the nozzle another jet is issued at M 0 c = 2.0. 
The freestream Reynolds number is 296,000. 

Figure 9 shows the inlet flow profiles of the ax- 
ial velocity, swirl velocity, radial velocity, pressure and 
density, which are taken from the experimental data 
of Ref. 12. The initial profiles are used as quasiaxi- 
symmetric profiles for the present computations. On the 
cylindrical boundary (side wall) of the flow at r = 3.5, 
freestream conditions are imposed corresponding to M x 
= 2.0. The initial conditions in the computational domain 
are also taken as those corresponding to the freestream 
conditions at M x = 2.0. The problem is solved using 
two types of exit boundary conditions at x = 7.0; first ex- 
trapolation of all five variables from the interior cell cen- 
ter and second using the Riemnn-invariant-type boundary 
conditions. 

II.l. Extrapolation from Interior Cell Centers 

Figure 10 shows snapshots of streamlines and Mach 
contours of the solution. The streamlines show multi- 
bubble breakdown at the early levels of time. These bub- 
bles develop due to the shock system formed at the nozzle 
exit in the vicinity of the centerline. It is noticed that a 
strong portion of the shock exists at the centerline which 
splits into two oblique shocks, one is a weak shock and 
the other is a strong shock. Behind the strong shock, the 
vortex breakdown bubbles exist. Thereafter, for t > 5, the 
oblique shocks move slowly in the downstream direction 
and breakdown the vortex bubble stays in its place. At t 
> 55, the shock system moves upstream and so does the 
vortex breakdown bubble. The slow motion of the shock 


system and the vortex breakdown bubble continues back 
and forth between these two locations. No vortex shed- 
ding has been captured during the computations of this 
case. It is also noticed that most of the exit points are 
continuously supersonic and hence no downstream effects 
exist with the exception of a very thin-layer around the 
centerline. 

II.2. Riemann Invariant Boundary Conditions 

Next, the boundary conditions at the exit are replaced 
by using the Riemann-invariant-type boundary conditions 
with pt = poe at the subsonic points. Figure 11 shows 
snapshots of the streamlines and Mach contours of the 
solution. By comparing the present solution with the 
previous case of Fig. 10, we see that there is very slight 
effect of the present boundary condition on the solution. 
This is understood since the subsonic region at the exit is 
very small and moreover, the exit boundary is far from 
the nozzle exit 

Figure 12 shows a blow-up of the Mach contours at t 
= 55 for the flow case of Fig. 10. The shock system near 
the nozzle exit is clearly seen. 

Concluding Remarks 

The unsteady, compressible full Navier-Stokes equa- 
tions are used to study the critical effects of the down- 
stream boundary conditions on the supersonic vortex 
breakdown. In the present study, two supersonic swirling 
flow cases are considered. The first one is for a super- 
sonic swirling flow in a configured circular duct, where 
four types of exit boundary conditions are used. The sec- 
ond one is for a supersonic swirling jet issued from a 
nozzle into another supersonic jet of lower Mach num- 
ber. This flow case is a computational simulation to the 
experimental case of Ref. 12. 

In the first flow case, we have shown that the 
downstream exit boundary conditions produce substantial 
changes in the vortex breakdown modes, particularly at 
advanced time levels. This is due to the effect of the sub- 
sonic downstream points on the flow upstream. In the sec- 
ond flow case, we have shown that the downstream exit 
boundary conditions do not produce substantial changes 
in the vortex breakdown modes. This is due to the very 
thin subsonic region at the exit and the large distance of 
the downstream exit from the nozzle exit 
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Figure 3. Blow-up of strea 
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Figure 7. Streamlines and Mach contours for a swirling flow 
with quasi-steady multibubble vortex breakdown, down- 
stream disk of r = 0.333. 





Figure 9. Quasi-axisymmetric flow profiles at x = 0.0 for super- 
sonic swirling jet from a nozzle. 














Figure 11. Streamlines and Mach contours for super- 
sonic swirling jet from a nozzle with low 
frequency almost single bubble vortex break- 
down, Riemann invariant B.C. 
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Figure 12. Blow-up of Mach contours at t = 55 for supersonic 
swirling jet from a nozzle with low frequency almost 
single bubble vortex breakdown, Riemann invariant 
B.C. 
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ABSTRACT- The unsteady, compressible, thin-hiyer Nav J“^“^ S ^ U ^! cy linder configura- 
SS—* •*— *• ’T^'lSriSr splitLg. finite-volume 
rinns The equations are solved using an p « s around a 5 ° semi-apex angle cone 

scheme. The computational ^ationsaiv" “^yu n drical afterbody of various length is added 
of unit length at various Reynolds numbe U y afterbody on the flow asymme- 

c ‘i-t 
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attention by researchers in the ^ in1p0II ance to dte dynamic stabdtty 

experimental fluid dynamics area ( • aircraft When flow asymmetry develops, it P uaas 

and controllability of missiles and fighter a recr aft- yawi „g, ^g and pitchmg moments that 

side forces, asymmetric lifting forces an system of the vehicle. 

might be larger than those avadable by the mmol y ady , thin-layer, compressible 

8 ln several recent papers by the present tmdtQre U.21. th agaric vonex flows 

Navier-Stokes equadons have been used ^™ulares^y ^ shapes . The emphasts of 

including their passive control, around con parameters which influence the asymmetric 

these papers was extensive computauonal stuped cost associated with dte solunon 
flow phenomenon and its passive control. Since * P rcsolution is vcry expensive, all the 

of ttaee-dimensional-flow problems with reastmab cal flow assumption. Such an assumption 

mTutauonal solutions were are in close proximity of eac 

reduces the problem solunon to that onoyo ^comcal p^ ^ magnitudc . Moreover, such 

^utioM^till provide\^ensive^ntos^^g^|^ e ^®^^J^ s ^^ l ^^[^^ I ^^ I ^ t ]| S ymmetric 

rX " Ic^Sudon near the solid .u^ry, 

2. FORMULATION AND C ® MPUT ^™s S A t ^n^y^om^sible, thin-layer Navier-Stokes 
The conservative form of the dunensmnles^ ,v* m used . The implicit, upwind, 

equations in terms of time-independent, baiy unsteady, compressible, thin-layer 

-iS tt m usci F “ 
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Fig ® 1. Aaymmetric flow solution around a con. of uni. length, short-duradon aid. slip. 
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Figure 3. Unsteady asymmetric locally-conical flow solution at different time steps within one 
cycle (cylinder axis is a time axis), At Iv . 


conditions a, specified a, die inflow fi — - 

is used at die outflow boundanes. The »ed « enforced and die 

solution. On the solid boundary, e no s ip Oie adiabatic boundary condition 

— ST- 

wdh L no- slip and no-peneuation conditions on die sobd boundary. 

3. COMPUTATIONAL APPLICATIONS AND DISCUSSIONS 

Circular Cone , enet h i s the characteristic length) 

■ A 5 rf ms ^5rS STi*. « considered by .he auihors in Ref. 1 for 

normal and axial directions, respecuvely, is gene y direction of the body using a 

a, axial stations. The grid is clusiemd algebnucaUy ” l(T> a, the axial suttion 

geometric series with minimum gnd s P acl " g ° „ d slze which was used for the 

of unit length. The cross-flow gnd size of 161x81 is the same gna 

locally-conical flow solutions of Ref. 1. =l8andRe=10 \ which are the same conditions 

With the flow conditions set at a - 20 , M«* ^ solution produces a symmetric 

as those of the locally-conical flow of Ref. 1, the asymmetric steady flow. Next, the 

steady flow, unlike the locally-comc so unor three-dimensional cone flow. In 
search is directed at obtaining asymmetric flow solutions tor tne ^ at 0 , 4 0°, M m 

Fig l 4 1 ’ Id an asymmetric and is nearly self-simil^over 

I « tie -e length. This solution is obtained using a short-durauon side-slip 
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Figure 4. 
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Figure 5. Asymmetric flow solution around a cone-cylinder configuratio 
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cone-cylinder configuration 0.5.10, a 
with experimental data (Ref. 5). 
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Circular Cone-Cylinder Configurations afterbody length on the flow asymmetry a 

To address the ^ ^ 

cylindrical afterbody of different . s solve d with the flow conditions of a = 40° Moo 

around the resulting cone-cylinder configura of ^ isolated unit-length cone of Fig. 1. 

- 1 4 and Re = 4x 10 6 , which are the same flow conch The of flow disturbance 

tophs of the cylindrical afrerberiy are choren “i e^Under configmanon of 1:1 (cone 
Tthe same short-duration 2°-side-shp disturbance^ on the. cone, in comparison 

length: cylinder length). Fig. 5 shows a very sno g ym ^ cylindrical afterbody as weU. It 
with the flow asymmetry of the isolated cone of Fig. h yli der configuration, subsonic 

^ould be noted ffiat inside the conical shock has an upstream effect, 

flow regions exist and hence the downstream fl ow asymmetry; the first is due to 

The cylindrical afterbody has dual effects which of ^ local angle of attack of 

dlT cone-cylinder juncture and the second is due toffie rn^ ^ ^ ^ ^ fl0 w , 

the leeward side of the cylinder. BoAaftoe e the asymmetry is strong and 

asymmetry. For the cone-cylinder configurations ot 

thHow becomes unsteady [Ref. 3]. ^ available experimental data. For this 

Next we show a comparison of the computed experimentally tested by 

purpo e, we consider the cone-cyhnder confi^tion of h L6 and the Reynolds 

Sum 5 . The configuration angle of C ytoder) is 6.6x10*. The cone semi- 
number based on the total configuration length ( ^ ^ 161x81x6 5. Figure 6 shows the 

apex angle is 9.5°. The problem is solved usuigag" total-pressure-loss contours and the 

surface-pressure coefficient along with the expe The computed and measured 

^ M^-number contours at the ax.ri ^o f O.475^ ^ ^ changes 

surface-pressure coefficient are in good agre tf axial stations at 0.475 and 0. . 

-r ^ ^ 

number for flows around a cone we have R lds flows. the spadal asymmemc 

changes sides in the downstream daemon. For th 7^ temporal asymmetnc flow 

flow Ivllops in a wavy manner, which ts Symmetry develops in a penodre 

development of the locaUy-comcal so^nons^ forcb ody, the flow asymmetry becomes 

manner By adding a cylindrical afterbody to a, the length of the cylindrical afterbody is 

stronger ^comparison wifl, ptally, dte compered resulri retd 

increased, the flow asymmetry becomes stronger an 
grid used are conclusively validated. 

1-994. The computations have been earned out using tn 
Research Center. 

REFERENCES r H ( 1990 ) “Prediction of Steady and Unsteady 

[1] Kandil, O. A., Wong, T-C. and ^; ^, aI ^a 90-0598. Also to appear in AIAA Journal, 
1 Asymmetric Vortical Bow Around Cones, AIAA 

Vol. 29, No. 9, 1991. 


1370 




j t • r w n99m “Numerical Simulation of Steady and 
Kandil, O. A., Wong, T-C. ^ £ SM Symposium on Non-Steady Fluid Dynamics, 

Unsteady Asymmetric Vortical Flows. ASItffi ympo Joimia i „f Fluids and Strictures. 

FED-Vol. 92, Toronto, Canada, pp. 1U0 - 1U 

1991. _ . » • C H Q991) “Three-Dimensional Navier-Stokes 

A^mertc ^oZlTtJ Con«-C^lmder ^nftgui^(^/’ ^ 9 ^ emcs 0 ver Slender 

ttttSZg from to d,3 for 

Landrum, E. J„ (1977) •Wmd-Ttatnel Prasum IJM NASA Technical 

a Series of Bodies of Revoludon at Angles of Attach trout 

Memorandum, NASA ^M-X-3558. Cylindcr wtth Lam.nar, Transitional, 

1 ^ «■ *>• no - "■ *• i492 - 14 "- 


1371 




Prediction of Steady and Unsteady 
Asymmetric Vortical Flows Around 
Circular Cones 

0. A. Kandil, T.-C. Wong, C. H. Liu 


Reprinted from 

AIM Journal 

Volume 29, Number 12, December 1991, Pages 2169-2178 

44MA 

A publication of the 

American Institute of Aeronautics and Astronautics. Inc. 
The Aerospace Center, 370 L Enfant Promenade, SW 
Washington. DC 20024-2518 




VOL. 29, NO. 12, DECEMBER 1991 


AIAA JOURNAL 


PRECEDING PAGE BLANK NOT FILMED 

2169 


Prediction of Steady and Unsteady Asymmetric 
Vortical Flows Around Circular Cones 

Osama A. Kandil* and Tin-Chee Wongt 
Old Dominion University , Norfolk , Virginia 23529 
and 

C. H. LiuJ 

NASA Langley Research Center , Hampton , Virginia 23665 


Steady and unsteady, supersonic asymmetric vortical flows and their passive control around circular cones are 
considered in this paper. These problems are formulated by using the unsteady, compressible, single and double, 
thin-layer, Navier-Stokes equations. The equations are solved by using an implicit, upwind, flux-difference 
splitting, finite-volume scheme, either in a pseudotime stepping or in an accurate lime stepping. An implicit, 
approximately factored, central-difference, finite-volume scheme has also been used to validate some applica- 
tions of the upwind scheme. Local conical flows are assumed for the computational applications presented in 
this paper. Steady asymmetric vortical flows have been predicted by using random and controlled disturbances. 
Unsteady asymmetric vortex-shedding flows have also been predicted, for the first time, using time-accurate 
solutions with two different computational schemes. Control of flow asymmetry has been demonstrated 
computationally by inserting a vertical fin in the leeward plane of geometric symmetry. 


Introduction 

I N the high angle of attack (AOA) range, the separated 
vortical flow from forebodies of missiles and fighter air- 
craft may become asymmetric, producing large abrupt 
changes in force and moment coefficients. These abrupt 
changes may exceed the available controllability and lead to 
missile and aircraft spin. Experimental studies of several re- 
searchers 1 " 11 have identified four distinct flow patterns about 
slender bodies through a wide AOA range and zero-degree 
side slip. The first pattern develops in the very small AOA 
range, where the flow is attached and the axial flow is domi- 
nant. In the intermediate AOA range, the crossflow becomes 
of the same order of magnitude as that of the axial flow, the 
flow separates on the leeward side, and a symmetric vortex 
pair is formed. As the AOA reaches a high range, the symmet- 
ric vortex pair becomes asymmetric, and the flows stay steady. 
For this asymmetric vortex-flow pattern to occur, it is not a 
necessary condition to have asymmetric separation lines on the 
leeward side of the body. The fourth flow pattern develops at 
a very high AOA range, where asymmetric time-dependent 
vortex shedding occurs either randomly or periodically, simi- 
lar to the von Karman vortex street in two-dimensional flows 
around cylinders. 

For isolated pointed forebodies, the onset of vortical flow 
asymmetry occurs when the relative incidence (ratio of AOA 
to nose semi-apex angle) exceeds a certain value; e.g., for a 
pointed circular cone, the relative incidence must be higher 
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than two. However, the relative incidence value is not the only 
determinable parameter for the onset of vortical flow asym- 
metry. The onset of vortical flow asymmetry is also a function 
of the freestream Mach number and Reynolds number and the 
shape of the body cross-sectional area as well. Asymmetric 
vortical flow and vortex shedding have also been documented 
for delta wings 12,13 at very high relative incidences and low 
subsonic regimes. 

For the critical values of the relative incidence, Mach num- 
ber and Reynolds number, and the shape of cross-sectional 
area, the symmetric flow is unstable. Any small flow distur- 
bance in the form of a transient side slip, acoustic disturbance, 
or similar source of disturbance causes flow instability that 
produces, depending on the flow conditions, either a steady 
asymmetric vortical flow or an unsteady asymmetric flow with 
vortex shedding. In this paper, we present an extensive compu- 
tational study of the steady asymmetric vortical flow and 
unsteady asymmetric flow with vortex shedding to address 
some of the influential parameters as the relative incidence 
and Mach number. 

As the experimental work shows, the mechanisms that lead 
to asymmetric vortex wake are not well understood. However, 
two mechanisms have been established for explaining the evo- 
lution of flow asymmetry. 5,6,9,10 The first mechanism applies 
to both laminar and fully turbulent flows. It suggests that flow 
asymmetry occurs due to instability of the velocity profiles in 
the vicinity of the enclosing saddle point that exists in the 
crossflow planes above the body primary vortices. 2,710 The 
second mechanism suggests that flow asymmetry occurs due to 
asymmetric transition of the boundary-layer flow either at the 
nose in the axial direction or on both sides of the body in the 
crossflow planes. For pointed slender bodies, the first mecha- 
nism produces higher side forces than those produced by the 
second mechanism. These results have conclusively been 
shown through the experimental work of Lamont 8,9 on 2-diam 
and 3.5-diam tangent ogive noses with cylindrical afterbody. 
An extensive review of the steady and unsteady vortex-induced 
asymmetric loads is given by Ericsson and Reding in Ref. 11. 

Several attempts have been carried out to computationally 
simulate asymmetric vortical flows around slender bodies of 
revolution. Early computational work on conical flows has 
been published in Refs. 14 and 15. Graham and Hankey 16 
presented the first three-dimensional Navier-Stokes computa- 
tions for asymmetric flow around a cone-cvlinder body at 30- 
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, . „ f attack i 6 freestream Mach number and 

?Sco“™ g'Sof 26 x 30 StoSSS; 

induced by ihe bod, 

thal triggers an instability of Ihe saddle point k 

(firs. mechanism tor ^fied .o pro- 

metry^can Jf *5ned 

banCC T^ a' m e, that was blown from one side of 
was in the form of a small jet t ^ was turned off> 

the body near the nose. • , showed that the flow 

,he numerical soluuon nntonunM^shonJd ^ ^ 
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time stepping, believed that the asymmetry was 

stayed stable thereafter. It is behevea in as a 

triggered by the mac hme ; ro ^ ^ nowfjeld ln a later 
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circular cone >n the angle of atta k g diameter . The 
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Formulation 

The three-dimensional compressible Jr^he dimen- 

i^rSJrT- i" tons are given b, 


where 


Q -j- j \p,pU\,pUi,pU^P e X 


( 2 ) 


£ m * inviscid flux 

= (\/j)[d k rE«Y 

= v/mun,^™ + a.rp.p«^«+ 

+ a,{"p.(pe + P)y-l'- m=1 ’ 2 ’ 3 ,3) 

( £ v): * viscous and head-conduction flux in I 2 direction 
= (l/./)[0,3*£ 2 T tl , d*S 2 T*2, (M-Tta-ff*)!' < 4 > 

(£v ) 3 s viscous and head-conduction fiux in S 3 direction 

= (l/i)lO,d*i } T/M, 3*£ 3 T* 2 , S k ^T k J, 3*i 3 (UnTjr„-9l)l 

i/ m =a*r«* (6) 

The first element of the three momentum elements of Eq. (5) 
is given by 


3*£ 3 r*i 


Afos/l 

5- kT 


(^a,£ 3 + 


dU] 

*a£ 


where 


«<> = d t £ J a*£ 3 . 


i at3 


an* 

as 3 


(7) 


( 8 ) 


The second and third elements of the momentum elements are 

2SSiS3. , SS5?. S’eS k (5) is gtven 
by 

a ,;|(«. r ..-«.)-^7 [* B ' + 4 W >“ J + “ !l 

i a(a 2 )H 
+ (-y- 1)P, ai 3 Jj 


(9) 


where 


W= 3„£ 3 Un 


( 10 ) 


r- c t*\ in the case of double thin-layer, Navier-Stokes 
For Eq. ( ), given by equations similar to Eqs. 

equations, the elements ar g m by ^2 jhe double 

sive control of flow asym y dicular to t he cone 

length, Revnolds number is defined as Re - p»L„L / 

P is related ,o .be rural energy pe, urn. 
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script or subscript m is a free index. The range of k , n , s , and 
m is 1-3, and d k ■ (d/d Xl ). 

Boundary conditions are explicitly implemented. They in- 
clude inflow-outflow conditions and solid-boundary condi- 
tions. At the plane of geometric symmetry, periodic condi- 
tions are used for symmetric or asymmetric flow applications 
on the whole computational domain (right and left domains). 
At the far- field inflow boundaries, freestream conditions are 
specified since we are dealing with supersonic Hows, whereas 
at the far- field outflow boundaries, first-order extrapolation 
from the interior points is used. On the solid boundary, the 
no-slip and no-penetration conditions are enforced; u , - u 2 
M, = 0 and the normal pressure gradient is set equal to zero. 
For the temperature, the adiabatic boundary condition is en- 
forced on the solid boundary. The initial conditions corre- 
spond to the uniform flow with u, = «i = «3 = 0 on lhe so11 

boundary. , r , 

For the passive control applications using a vertical tin in 
the leeward plane of geometric symmetry, solid-boundary 
conditions are enforced on both sides of the fin. 

Highlights of Computational Schemes 

The first computational scheme used to solve the unsteady 
compressible, single or double thin-layer, Navier-Stokes equa- 


tions is based on the Roe inviscid flux-difference splitting 
scheme In this scheme, the Jacobian matrices of the inviscid 
fluxes. A, = (dE s /dq), s = 1-3, are split into left and right 
fluxes according to the signs of the eigenvalues of the inviscid 
Jacobian matrices. Flux limiters are used to dampen the nu- 
merical oscillations in regions of large changes of the gradients 
of the flowfield vector. The viscous and heat transfer terms 
are centrally differenced. The resulting equation is solved by 
using approximate factorization in the {', £ , and ? direc- 
tions. The computational scheme is coded in the computer 

program CFL3D. . 

The second computational scheme is an implicit, approxi- 
mately factored, centrally differenced, finite-volume scheme." 
Added second-order and fourth-order dissipation terms are 
used in the difference equation on its right-hand side terms, 
which represent the explicit part of the scheme. The Jacobian 
matrices of the implicit operator on the left-hand side of the 
difference equation are centrally differenced in space, an 
implicit second-order dissipation terms are added for the 
scheme stability. The left-hand side operator is approximately 
factored, and the difference equation is solved in three sweeps 
in the £‘, e, and £ } directions, respectively. computa- 
tional scheme is coded in the computer program 1CF3D. the 
ICF3D code is used to verify some of the applications ot the 





Fig. 1 


Symmetric flow solutions for a circular cone, a 


10 deg, M. = 1.8, Re = 10 s (validation case). 



Symmetric 


— CFL3D 




due to random disturbances, a — 20 deg, 


= 1.8, R, = 10 s (validation case). 


Fig. 2 Steady asymmetric flow solutions for a circular cone 
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CFL3D code; namely the cases of Figs. 1 and 2^ For the 
problem of passive control of flow asymmetry, the double 
thin-layer, Navier-Stokes equations have been solved using the 

C *Since ^applications in this paper cover local-conical flows 
only the three-dimensional scheme is used to solve for locally 
conical flows. This is achieved by forcing the conserved com- 
ponents of the flow vector field to be equal at two planes of 
5 =0.95 and 1.0. The validity of local-conical-now assump- 
tion is discussed in the next section. 

Validity of the Local-Conical-Flow Assumption 

The solutions presented in this paper are called local-conical 
solutions which are obtained by equating the conserved com 
ponints of the nowfield vector, in the three-dimensional 
scheme on two crossnow planes that are in close proximity to 
Sh mher at a selected location. Once this location is specified 
(" =10 in the present applications), the now Reynolds num- 
ber is determined and the time scale, for time-accurate solu- 
tions is also determined. The resulting solution is a local-ton- 
ical solution at the specified location. It is not a global-comca 
solution The locally conical equations can be shown by con- 
sidering the conservative form of the Navier-Stokes equations 
in the Cartesian system 

(13) 


dq d{£ Ev)i q 

Tt + d x, 


i = 1 - 3 


By introducing the conical coordinates 


V\ ■ 


*i 

aV 


Tj *> — — , XjXj 

* 


(14) 


and using the chain rule to express Eq. (13) in terms of the 
conical coordinates, we get 

^ + A ( £_£„) i+ -^(E-£ v ): 

6 dt drfi ^2 


+ 52 — (£ - £,.), + 2(7- ly) = 0 
e 2 dm 


(15) 


where 


8 = V 1 +»); + »): 

£,=£,- t),£ 3 , £; = £: _ ViE} 

£, = £j + i)|£i + 

7 = £j 

£ vl = £, 1 - i?i £.3 
£ l2 = £,; - I)2£v3 
£,. 3 = £,3 + 7 i £ >i + rj 2 £v2 
/, = £,3 


(16) 


(17) 


The conical How condition requires that the flow variables be 
independent of the coordinate *. If this condition is imposed 
in Eq. (15), by dropping the derivatives with respect to tj,, the 
equation reduces to 

0! d J. + JL ( £_£ V ), + i-(£-£.h + 2 (/-A) = 0 (18) 

8 dt 3t), drjj 

It is clearly seen that Eq. (18) still has rj, dependence in the 
unsteady term and the viscous and heat-nux terms (one can se 
the explicit dependence of the viscous and heat-flux terms on 
' bv transformmg the elements of these vectors to the conical 


coordinates). Hence, Eq. (18) is not self-similar, and there- 
fore it does not represent a global-conical flow. Howwer. m 
is set equal to a constant e. then one can consider Eq. (18) to 
represent a local-conical flow around ij, -c. The result! g 
solution using Eq. (18) with tj, = c represents a local-conical 
solution with a Reynolds number and a time that are scale > 
the constant c. It should be noted that if the flow is steady and 
: S. Eq. (18) becomes self-similar, and bene. rep- 

resents a global-conical How. In the present paper, indi- 
rectly solve Eq. (18) at a fixed location of un.ty^ This 
achieved in the three-dimensional flow equation, Eq. (1), 
equating the elements of the flowfield vector at two planes 
close proximity to each other. In this paper, we selected these 
be located a, a - 0.95 and 1.0 In otto nqme, c.1 
experiments, we use the plane locations at x = 0.995 and 1 .0 
The results of these experiments were in excellent agreement 
with those of the present paper. 

Computational Studies 

Supersonic nows about a 5 deg semiapex angle circular cone 
at a Reynolds number of 10 5 have been considered. A grid o 
161 x 81 points in the circumferential and normal directions s 
used throughout the present applications. The grid i g 
ated by using a modified Joukowski transformation w ith 
geometric series for the grid clustering ; neat r the ^ 
boundary. The minimum grid length is 10 at the 
boundary, and the maximum radius of the computation^ 
domain is 21r , where r is the radius of the circular cone at the 
axial station of unity. 

Steady Symmetric Flows 

Figure 1 shows steady symmetric vortical-flow solutions for 
the circular cone at 10 deg angle of attack and 1.8 freestream 
Mach number. In the figure, we show comparisons of the 
results of the CFL3D and 1CF3D codes. The results include 
the 11 residual error versus the number of iterations the cross- 
flow velocity, the total-pressure-loss contours and the sur 
face-pressure coefficients. It should be noted here that the 
angle 8 in the C p figure is measured from the leeward plane o 
geometric symmetry in the clockwise direction^ axemen 
of the results of the two code is excellent and the results are 
in full agreement with those of Siclari and Marconi. 

Steady Asymmetric Flow 

Round-Off and Truncation Error Disturbances 

The cone angle of attack is increased to 20 deg while all the 
other How conditions are kept fixed. Figure 2 shows the result 
of the CFL3D and 1CF3D codes. In the residual error figure, 
the CFL3D code shows that the residual error drops 10 orders 
of magnitude within 2500 iteration steps. Thereafter, the error 
increases by six orders of magnitude. The How is symmetric 
during this 5000 iteration steps. Next, the error drops do y 
another six orders of magnitude and stays constan f ° r 
iteration steps. The now becomes asymmetric and stable. The 
1CF3D code shows that the residual error drops ive or ers 
magnitude in the first 3000 iteration steps, increases two or- 
ders of magnitude in the next 2000 iteration steps, and then 
drops down by three orders of magnitude within the next 5000 
iterations. The flow solution goes through a symmetric un- 
stable solution and then to an asymmetric stable soluuon. T ^ e 
pressure-coefficient figure for the two ^es ,s .the same over 
the full range of the circumferential angle 6. The suctio 
nressure in the range of 8 = 0-90 deg is lower than that of the 
range of 8= 270-360 deg. The crossflow velocity and tota - 
pressure-loss contours for the two codes are also in excellent 
agreement They show the nature of the How asymmetry an 
its details. The results are in complete agreement with those of 

RC Since the residual error of the CFL3D code is much smaller 
,h,r,b« qf .he ICF3D code after .he firs, 
disturbance Iba, .r.gsered the asymmetry m .be fra : «* ' » 
attributed to the machine round-off error, while the di 
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Fig. 3 



Steady asymmetric flow solutions for a 


circular cone due to ±2 deg transient side slip. 


a = 20 deg* Af« = 


R e = 10 5 . 


Fig. 4 



Steady asymmetric flow solutions for a circular cone due to 


± 5 deg transient side slip, or = 20 deg. Mm 


= 1.8. Re = 10 s - 


bance that triggered the asymmetry in the second code is 
attributed to the truncation error of the scheme (since there is 
a bias due to the spatial marching direction). Both distur- 
bances are random in nature. However, irrespective of the 
source of disturbance, the final asymmetric stable solution is 
the same. 

Controlled Transient Side-Slip Disturbances 

In Figs. 3 and 4, we show steady asymmetric flow solutions 
due to transient side-slip disturbances of ± 2 and ± 0 5 deg. 
The residual-error figures show a drop of seven orders of 
magnitude in the first 2000 iterations. At this step, a side-slip 
disturbance is imposed for six iteration steps, then it is re- 
moved. Irrespective of the magnitude or the sign of the side- 
slip disturbance, the residual error increases by six orders ol 
magnitude, then it drops down very rapidly. A stable asym- 
metric flow solution is obtained. The asymmetric solutions 
corresponding to the ± 2 deg side-slip disturbances are mirror 
images of each other, as can be seen from the figures of the 
surface-pressure coefficient, crossnow velocity, and total- 
pressure-loss contours. The corresponding asymmetric solu- 
tions with the ± 0.5 deg side-slip disturbances are exactly the 
same as those of the ± 2 deg side-slip disturbances. Moreover, 


the final asymmetric solutions of the ± 2 deg and ± 0.5 deg 
side-slip disturbances are the same as those of Fig. 2. 

Again, this numerical experiment shows that the same phys- 
ical How asymmetry is obtained. 


Insteady Asymmetric Vortex Shedding 

In the present case, the angle of attack is increased to 30 deg 
nd all the other now conditions are kept the same as those ot 
he cases above. Figure 5 shows the results of this case. 

Here we show the history of the residual error and the lift 
oefficient up to the 15,700 time step. First, pseudo-time 
tepping was used up to 10,000 iterations, and the solution was 
nonitored every 500 iterations. The solution showed that the 
isymmetry was changing from the left side to the right side, 
vhich indicated a possibility of unsteady asymmetric vortex 
shedding. The residual error was also oscillating. The compu- 
lations were repeated starting from the 3,500 iteration step 
jsing time-accurate calculations with At - 10 . The resM *u*'- 

error and lift-coefficient figures show the time history of the 
solution. It is seen that the residual error and the lift coeffi- 
cient show a transient response that is followed by a periodic 
response. Figure 5 shows also snapshots of the time history or 
the solution for the total-pressure-loss contours and surface- 
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limiters are turned on. The resulting numerical dissipation in 
the FVS is large enough to dampen the random disturbances 
of the flow solution. By turning off the flux limiters in the 
FVS scheme, the random disturbances can grow, producing 
the asymmetric unsteady vortex shedding. This also shows 
that the FDS scheme, even with the flux limiters turned on, is 
less dissipative than the FVS scheme. These results conclu- 
sively explain the erroneous claim of steady flow made by the 
previously mentioned researcher. 

Steady Asymmetric How at Different Mach Numbers (Effect of M„) 
Figure 7 shows the effect of the freestream Mach number 
= 2.2, 2.6, and 3.0) on the convergence history, surface 
pressure! crossflow velocity, and total-pressure-loss contours 
for the circular cone at 20 deg angle of attack. At = 2.2, 
the residual error shows that the stable asymmetric flow is 
obtained within the same number of iterations as that of the 
A/„ = 1.8 case. At = 2.6, the residual error shows that the 
stable asymmetric flow is obtained after a large number of 
iterations. And at A/. = 3.0, no asymmetric How was cap- 
tured, the flow stayed symmetrically stable. The surface pres- 



sure figures show that the asymmetry gets weaker as the Mach 
number is increased. This conclusion is clearly seen from the 
crossflow velocity and the total-pressure-loss figures. It should 
be noted that since the nature of disturbance is random. How 
asymmetry changes sides as the Mach number increases until it 
disappears. 


Passive Control of Flow Asymmetry 

Figure 8 shows the passive control of flow asymmetry by 
inserting a vertical fin in the leeward plane of geometric sym- 
metry. The fin height is equal to the cone local radius r. Here, 
the double thin-layer, Navier-Stokes equations are used to 
obtain these results. The flow Mach number is kept at 1 .8 and 
the angle of attack is 20 deg. The flow is completely symmetric 
as can be seen from the figures of the surface-pressure coeffi- 
cient, total-pressure-loss contours, and crossflow velocity. A 
blow-up of the cross-flow velocity at the fin-cone juncture 
shows two corner recirculating bubbles of exactly the same 
size. This case has been obtained after 24,000 iteration steps. 
Again, this is the first time such a computational simulation of 




Fig. 6 Unsteady asymmetric flow solution (using flux-vector splitting) with vortex shedding for a circular cone during 
a = 30 deg. = 1.8, R, = 10 5 , At = 10 3 . 


periodic flow response; 
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Fig. 8 Control of asymmetric flow of a circular cone using a vertical fin in the leeward geometric plane of symmetry, a = 20 deg, - 1.8, 
Re = 10 s , h = r (r m radius of circular section). 


the passive control of the flow asymmetry has been presented. 
The results are in full agreement with Stahl’s experimental 
study. 21 

Concluding Remarks 

This paper presents extensive computational study and sim- 
ulation of steady and unsteady asymmetric vortex flow around 
circular cones. A systematic study has been carried out to 
show the effects of angle of attack and Mach number. The 
study shows that the flow asymmetry is independent of the 
type or level of the disturbance. For the controlled transient 
side-slip disturbance, the solution is unique. For the uncon- 
trolled random disturbance, the solution is also unique with 
the exception of having the same asymmetry changing sides on 
the cone. It conclusively shows that periodic vortex shedding 
has been captured at larger angles of attack. The unsteady 
asymmetric vortex-shedding solution has been substantiated 
by using two different computational schemes. It also shows 
that as the Mach number increases, the vortex flow asymmetry 
gets weaker until it disappears. The possibility of passive 
control of flow asymmetry has also been demonstrated. Many 
of the cases presented here are obtained for the first time, in 
particular, the asymmetric vortex shedding cases and the cases 
of passive control of flow asymmetry. 
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The unsteady compressible, thin-layer, Navier-Stokes (NS) equations are solved 

i ♦ and unsteady asymmetric, vortical laminar flow around cones at high 

simulate steady a^unsteimy^ asym^^ ^ tlons are OTlved by using an implicit. 

,,d iffereS splitting (FDS), finite-volume scheme. The locally con.cal flow 
upwind, flux d P } obtained by forcing the conserved components of 

SE vSi?!b?S^I a. 1*0 axial stations located a. 0-95 and 10. ComP“mio„> 

examples cover scfiufion ar<m^d t^e drcularcone 

control using side strakes. The u M vector solitting (FVS) scheme with the 

round- and sharp-edged delta wings. 


1. INTRODUCTION 

At high angles of attack, flow separations from the forebodies of missiles and fighter 
aircraft may become asymmetric resulting in side forces, yawing moments and ro ing 
moments which are in many instances, sufficiently large to trigger missile and aircraft 

^tn'^xperimenta^ studies have shown that it is not necessary for the separation lines 
spin, t-xperimemai mu e Pnar ated flow to be asymmetric (Kenner & 

Chapmrw^'pea'k art. 1979; Lamom 1980, 1982). These studies have also shown 
AafuTeady asymmetric flow with vortex shedding may be either random or penodic 
where the latter ri similar to the Kfirmfin vortex street m two-dimensional flows aroun 

"“Thronse. of flow asymmetry occurs when the relative incidence (ratio of angle of 
attack to semi-apex angle) of pointed forebodies exceeds certain critical values At the 
critol vatos of relative incidence, flow asymmetry develops due to natural and/or 
forced disturbances The origin of natural disturbances may be a transient side slip, a 

Kna ffi'oneoni ihedeSSbte parameters for the onset of flow 
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have also been documented for sharp-edged delta wings at very high relative incidences 
(Shanks 1963; Rediniotis et al. 1988). 

The mechanisms which lead to steady and unsteady asymmetric vortical flows past 
wings and bodies at high angles of attack and zero side slip are not well understood. 
The experimental studies of these phenomena by several investigators (e.g. Keener & 
Chapman 1977; Lamont 1982; Skow & Peake 1982; Peake & Tobak 1982) propose two 
mechanisms for explaining the origin of flow asymmetry. The first mechanism suggests 
that the asymmetry occurs due to instability of the velocity profiles in the vicinity of the 
enclosing saddle point which exists in the cross-flow planes above the body primary 
vortices. The second mechanism suggests that the asymmetry occurs due to asymmetric 
transition of the boundary-layer flow at the apex, either in the axial direction or on 
both sides of the body in the cross-flow plane. 

Very recently, several attempts have been carried out to computationally simulate 
asymmetric vortical flows around slender bodies of revolution. In a paper by Marconi 
(1988), the Euler equations are used along with a “forced separation model”, which is 
introduced by Fiddes (1989), to solve for supersonic flow past a circular cone. The 
pseudo-time stepping is carried out until the residual error reaches machine zero while 
the flow is symmetric. Proceeding with the time stepping, vortex-flow asymmetry is 
obtained and stays stable thereafter. It is believed that the asymmetry is triggered by 
the machine round-off error, which acts as a disturbance to the saddle point in the 
flowfield. This work shows that the first mechanism of asymmetric vortex flow is 
basically an inviscid mechanism. In a recent paper by Siclari & Marconi (1989), the full 
Navier-Stokes equations are used to solve for supersonic asymmetric flows around a 
5°-semi-apex angle cone over a wide range of angles of attack. 

Kandil et al. (1990a) used the unsteady, thin-layer Navier-Stokes equations along 
with two different implicit schemes to simulate asymmetric vortex flows around cones 
with different cross-sectional shapes. The numerical investigation focuses on a 
5°-semi-apex angle circular cone and locally conical flow is assumed. The first scheme is 
an implicit, upwind, flux-difference splitting, finite-volume scheme and the second one 
is an implicit, approximately factored central-difference, finite-volume scheme. Keep- 
ing the Mach number and Reynolds number constant at 1-8 and 10 5 , respectively, the 
angle of attack is varied from 10 to 30°. At a = 10°, a steady symmetric solution is 
obtained and the results of the two schemes are in excellent agreement. At a = 20° and 
irrespective of the type or level of the disturbance, a unique steady asymmetric solution 
is obtained and the results of the two schemes are in excellent agreement. Two types of 
flow disturbances are used: a random round-off error or a random truncation-error 
disturbance, and a controlled transient side-slip disturbance with short duration. For 
the controlled transient side-slip disturbance the solution is unique, and for the 
uncontrolled random disturbance the solution is also unique with the exception of 
having the same asymmetry changing sides on the cone. At a = 30°, an unsteady 
asymmetric solution with vortex shedding is obtained, and the vortex shedding is 
perfectly periodic. Next, the angle of attack is kept fixed at 20° and the Mach number is 
increased from 1-8 to 3 0 with a step of 0-4. The solutions show that the asymmetry 
becomes weaker as the Mach number is increased. The flow recovers its symmetry 
when the Mach number reaches 3*0. Passive control of the flow asymmetry has also 
been tentatively demonstrated by using a fin on the leeward side of the body along the 
plane of geometric symmetry. 

Experimental research efforts have also been directed to control asymmetric flows 
for eliminating or attenuating the asymmetric forces and the resulting moments by 
using either passive-control or active-control methods. Passive-control methods include 
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the use of a vertical fin on the leeward side along the ^**^*2?* 
(Stahl 1989), the use of fixed or mw-able cr oss-section, from a 
1989), or the use of a rotatable '"7* tip (Moskovirz e, at. 1990). 

circular shape at its base to an elliptic s p ith V{ir ; om hlowinp 

Active-control methods primarily ^“turfa'c^Ng 1 W). Computational simulations 
r £ £££ iVX thOeffecriveOesTofOmth pais ive (Kandi, e, nf. 1990a, and 

active control methods (Travella et al. l" 0 )* . . . Navier-Stnkes 

In this paper, the unsteady, compr^l ble am. finite . 

equations are used, along wrth rm ^rmpha P e(rjc vortex flows around 

volume scheme to solve for steady ^ and unst dy * dreolaI cone and its 

cones. The steady rcsults_.nc h.de asj nnmnmc^ etric vortex shed d,ng 

control using side strakes. The unsteady results mwuu J , , f the 

around circular elliptic- and diamond-section cones. The unsteady results for the 
circular cTneare verified by using the FVS with the thin-layer Navrer-Stokes equatrons 
and the FDS with the full Navier-Stokes equations. 


2. FORMULATION 


2.1. Governing Equations 

The three-dimensional compressible viscous flow around the body is governed by the 
conservative form of the dimensionless, unsteady, compressible, thin-layer Navier 
Stokes equations. In terms of time-independent body-conformed coordinates § , § 
and i 3 , the equations are 

B<2 9 E s 9(^)_ n (1) 

aT + 3r 3S 1 ’ 


where 

q = 1 = - [ p> pUl , pu 2 , pu 3 , pe]'; 

J J 

E m = inviscid flux = - [3*f m £*] < 

= j [pU m , pui U m + 9i % m p, pu 2 U m + 3 2 % m p, pUiU m 

+ 3 3 ^ m p. (P* + P)Vm\, m = 1, 2, 3; 

(£ u ) 3 = viscous and heat-conduction flux in § 3 direction 

= j [ 0 , 3*| 3 t*„ 3*£ 3 t* 2 , 3*! 3 t* 3 , 3 k | 3 (u„T*„ - <?*)]'; 

u m = 3*r , «*. 

The first element of the three momentum elements of equation (4) is given by 



( 2 ) 


(3) 

(4) 

(5) 


( 6 ) 
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tf> = 9*S 3 * 3*S 3 , 


V = i3*t 3 ||f. 


( 7 ) 


The second and third elements of the momentum elements are obtained by replacing 
the subscript 1, everywhere in equation (6), with 2 and 3, respectively. The last 
element of equation (4) is given by 


- <lk) '■ 


M®/x 

Re 




1 3 (a 2 ) 

23§ 3 “ J “ 5 (y-l)Pr 9£ 3 


]}• 


( 8 ) 


where 

W = 9 „| 3 u„. (9) 

The reference parameters for the dimensionless form of the equations are L, a®, L/a„, 
p. and Hoc for the length, velocity, time, density and molecular viscosity, respectively. 
The Reynolds number is defined as Re = and the pressure, p, is related to 

the total energy per unit mass, e, and density, p, by the gas equation 

P * (y ~ 1 )p(e ~ 2 u„u n ). (10) 

The viscosity, p, is calculated from the Sutherland law 

'‘"’""(f+l)’ C = 0-4317, (11) 


and the Prandtl number Pr = 0-72. 

In equations (1)-(10), the indicial notation is used for convenience. The subscripts k 
and n are summation indices, the superscript or subscript s is a summation index and 
the superscript or subscript m is a free index. The range of k, n, s and m is 1-3, and 
d k = d/dx k . In equations (1)— (11), ix„ is the Cartesian velocity component, U m the 
contravariant velocity component, r m the Cartesian component of the shear stress 
tensor, q k the Cartesian component of heat flux vector, a the local speed of sound and 
M„ the free-stream Mach number. 


2.2. Boundary and Initial Conditions 

Boundary conditions are explicitly implemented. They include inflow-outflow condi- 
tions and solid-boundary conditions. At the plane of geometric symmetry, periodic 
conditions are used for symmetric or asymmetric applications on the whole computa- 
tional domain (right and left domains). Since we are dealing with supersonic flows, at 
the far-field inflow boundaries, free-stream conditions are specified, and the conical 
shock is captured as part of the solution. At the far-field outflow boundaries first-order 
extrapolation from the interior points is used. On the solid boundary, the no-slip and 
no-penetration conditions are enforced, Ui = u 2 = u 3 = 0, and the normal pressure 
gradient is set equal to zero. For the temperature, the adiabatic boundary condition is 
enforced on the solid boundary. The initial conditions are set equal to the free-stream 
conditions with u 2 = u 2 = «3 = 0 on the solid boundary. 

For the passive control application using side strakes, solid-boundary conditions are 
enforced on both sides of the strake. 


3. COMPUTATIONAL SCHEMES 

The principal computational scheme used to solve the governing equations is an 

implicit, upwind, flux-difference splitting, finite-volume scheme. It employs the 
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flux-difference splitting scheme of Roe. The Jacobian matrices of the mviscid 
flux-differences are split into left and right flux differences according to the signs of the 
eigenvalues of the inviscid Jacobian matrices. The smooth flux limiters are used to 
eliminate oscillations in the shock region. The viscous and heat-flux terms are centrally 
differenced. The resulting difference equation is solved using approximate factorization 
along the f 1 , § 2 and § 3 directions, respectively. The scheme is third-order accurate in 
space and first-order accurate in time. The computational scheme is co e m t e 
computer program “CFL3D.” Details of the scheme are given by Rumsey & Anderson 

(1988). 

The second scheme is an implicit, approximately factored, central-difference 
finite-volume scheme. Added second-order and fourth-order dissipation terms are use 
in the difference equation on its right-hand side terms, which represent the explicit part 
of the scheme. The Jacobian matrices of the implicit operator on the left-hand si e o 
the difference equation are centrally differenced in space, and implicit second-order 
dissipation terms are added for the numerical stability. The left-hand side operator is 
approximately factored, and the difference equation is solved in three sweeps along the 
I 1 , § 2 and f 3 directions, respectively. The computational scheme is coded in the 
computer program “ICF3D”. Details of the scheme are given by Kandil & Chuang 
(1989). The ICF3D code is used to verify some cases which are solved by the CFL3D 
code. In this paper, the ICF3D code is used to validate the steady asymmetric flow 
case. 

A third computational scheme is used to validate the unsteady asymmetric vortex 
flow around the circular cone at a 30° angle of attack. This scheme is the flux-vector 
splitting scheme which is based on the van Leer flux-vector splitting scheme (Rumsey 
& Anderson 1988). This optional scheme is also coded in the computer program 
“CFL3D”. 

Since the applications in this paper cover locally conical flows only, the three- 
dimensional codes are used to solve for locally conical flows at the axial station X\ = 1. 
This is achieved by forcing the conserved components of the flowfield vector, q, to be 
equal at two planes located at jc a = 0-95 and 1-0. The concept of locally conical flow is 
explained in the next section. 

4. LOCALLY CONICAL FLOWS 

Locally conical solutions of the thin-layer or full Navier-Stokes equations are obtained 
using one of two methods. In the first method, the governing equations are 
transformed using the conical-coordinate transformation. Invoking the conical flow 
condition which requires that the flow variables be independent of the radial distance 
(or axial distance, depending on the transformation) from the cone apex, equating the 
radial distance (or axial distance) which appears in the transformed equations to a 
constant (equals to unity in most of the present locally conical solutions), the resulting 
equations are solved on one spherical (or cross-flow) surface. In the second method, 
the three-dimensional flow equations are solved on two spherical (or cross-flow) 
surfaces which are located in the very near proximity of a constant radial (or axial) 
distance. 

During the pseudo-time or accurate-time stepping, the flowfield vector is forced to 
be equal at the corresponding grid centers on the two surfaces. This method is used in 
the present paper to obtain locally conical solutions. The resulting solutions from these 
two methods are the same locally conical solutions. These solutions correspond to the 
specified radial (or axial) distance and hence they change as the radial (or axial) 
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distance is changed. The reason behind that is simply because the transformed 
equations, according to the first method, are not self-similar and hence they are not 
globally conical. This is shown below by developing the transformed equations of the 
first method. Considering the unsteady, compressible, Navier-Stokes equations in the 
Cartesian coordinates, 


Bq 9 (E - EJ, 

9 1 dx, 


introducing the conical coordinates, 


*1 *2 

1 h=-> r ?2 = ~, 

Xj x 3 


i = 1-3, 


T}1 = XiX h 


and using the chain rule, equation (12) is transformed to 

-^ + ^(E-E v ) l +^(E-E v ) 2 + ^^(E-^) 3 + 2(i-l v ) = 0, 
m 9 1 9 rj, 9 rj 2 m 9tj 3 

where 


m = Vl + rjT+ vl> E, = Ej - tj,E 3 , 

£2 = E — TJ2E3, £3 = E 3 + tj.E, + 7I2E2, 

i = £ 3 , £„, = E„, - T/jE., 3 , 

£t)2 = E„ — Tj2^v3’ £«3 = E„3 + r}]E„] + f)2 Et)2i 

f- = £v3- 


( 12 ) 


(13) 


(14) 


(15) 


The conical flow condition requires that the flow variables be independent of the 
coordinate rj 3 (radial distance). Invoking this condition in equation (14) by dropping 
the derivatives with respect to rj 3 , equation (14) reduces to 


-^ + ^-(E-t v ) l +^-(E-E v ) 2 + 2(\-l) = 0. (16) 

m 9 1 977 , 9ij 2 

It is obvious that the unsteady term includes tj 3 . Moreover, the viscous terms 
9 £ v1 /9t7 1 , dE v2 /dr ]2 and include q 3 , and hence equation (19) is not self-similar. The 
explicit dependence of the viscous terms on r\ 3 can be shown through one of the 
elements of these vectors. For example, we consider 


9 . Moc 9 

-(r xx -r] x x xz ) = — — 
9tji R 3 9rh 


M-m 9 


Re*?3 9»J 


;WG 


( 2 *±i 

2 

/9u, 

9 u 2 

l z 

l dU3 \ 

/ 9«i 

V 9x 

~ 3 

\ 9jc 

By 

Bz) 

\ oz 


9 a. 

L 

Bu, 

2 Bu 2 


") 

o> 

cs 

+ ^1 

Bij 2 

3 9»? 2 

3 9rj, 


2 9m 3 
rr/ 2 — 

3 drj 2 


)]■ 


(17) 


Thus, the unsteady term and viscous terms are scaled by the radial distance rj 3 and 
equation (16) does not represent a globally conical flow. The best to be done to make 
use of this equation is to select a constant value for ij 3 , and solve the resulting equation 
for what we call “locally conical flow”. If q 3 is assigned another constant value, the 
resulting equation will have another scale for the unsteady term and viscous terms. It is 
concluded that equation (16) becomes globally conical if the unsteady term and viscous 
terms vanish, and hence only the steady Euler equations are globally conical. 
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5.1. Steady Asymmetric Flow Around a Circular Cone 

A grid of 161 x 81 points in the circumferential and normal directions, respectively is 
used. The grid is generated by using a modified Joukowski transformation with a 
geometric series for the grid clustering near the solid boundary. The minimum grid 
spacing at the solid boundary in the normal direction is 10 -4 (the characteristic length 
is the conical station at x x = 1). A 5°-semi-apex angle cone at angle of attack, a, of 20°, 
free-stream Mach number, M., of 1-8 and Reynolds number. Re, of 10 s is used. The 
maximum radius of the computational domain is 21 r, where r is the cone cross-section 
radius at the axial station x, = 1. 

Figure 1 shows the residual error versus the number of iterations, the surface- 
pressure coefficient, the cross-flow velocity and the total-pressure-loss contours for the 
CFL3D and ICF3D codes. In the residual error figure, the CFL3D code shows that the 
residual error drops ten orders of magnitude, to machine zero, within 2,500 iteration 
steps and the solution is still symmetric. Thereafter, the error increases by six-orders of 
magnitude and slightly asymmetric solutions are obtained. The flow is symmetric 
during the first 5,000 iteration steps. Next, the error drops down by another six orders 
of magnitude, to machine zero again, and stays constant, and a stable steady 
asymmetric solution is obtained. It should be noted here that when the residual error 
first reaches machine zero, the machine-round-off error acts as a random asymmetric 
disturbance to the critically symmetric solution. Thereafter, the disturbance grows 
spatially, producing the asymmetric solution. The ICF3D code shows that the residual 
error drops five orders of magnitude in the first 3,000 iteration steps, increases two 
orders of magnitude in the next 2,000 iteration steps, and then drops down by three 
orders of magnitude within the next 5,000 iterations. The flow solution goes through a 





Figure 1. Comparison of steady asymmetric flow solutions around a circular cone, a = 20°, NL = 1-8, 

Re = 10 s . 
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symmetric unstable solution and then to the asymmetric stable solution. The pressure 
coefficient curves for the two codes are in excellent agreement over the circumferential 
angle 6, which is measured from the leeward plane of symmetry. The suction pressure 
in the range of 0 = 0-90° is lower than that in the range of 6 = 270-360°. The 
cross-flow velocity and total-pressure-loss contours for the two codes are in excellent 
agreement. They show the nature of the flow asymmetry and its details. 

Since the residual error of the CFL3D code is much smaller than that of the ICF3D 
code, the disturbance which triggered the asymmetry in the first code is attributed to 
the machine round-off error, while the disturbance which triggered the asymmetry in 
the second code is attributed to the truncation error of the scheme. Both disturbances 
are random in nature. However, irrespective of the source of disturbance, the final 
asymmetric stable solution is unique. Kandil et al. (1990a) have shown that the solution 
is still unique if another source of disturbance is applied for the same critical flow 
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conditions; e.g., a 2° or 0*5° short-duration side-slip disturbance produces the same 
flow asymmetry. 

In Figure 2, we show the effect of the minimum grid size on the flow asymmetry for 
the same flow conditions of the previous case. The figure shows the residual error 
versus the number of iterations, the total-pressure-loss contours and the surface- 
pressure coefficient for = 10~ 3 , 10 -4 and 1(T 5 at the cone surface. The histories 
of the residual errors are qualitatively of similar behavior. The total-pressure-loss 
contours show unique solutions with the exception of having the asymmetry changing 
sides. This is understood due to the random nature of the disturbance — a machine 
round-off error. The surface-pressure-coefficient curves also show unique solutions. 

With all the numerical experiments given above, it is conclusively proven that the 
asymmetric solution is not scheme-, numerics- or disturbance-dependent. 

5.2. Passive Control of Asymmetric Flow Using Side Strakes 

Figure 3 shows the results of passive control of flow asymmetry around the circular 
cone of Figure 1 by using side strakes of height equal to 0*3 r. The iteration histories of 
the residual error, lift coefficient and side-force coefficient show the attenuation of the 
flow asymmetry and the final stable symmetric solution. The surface-pressure- 
coefficient curve, the cross-flow velocity and the total-pressure-losses contours show the 
final symmetric solution and the symmetric vortices associated with this controlled flow. 




Figure 3. Passive control of asymmetric flow around a circular cone using strakes, a = 20°, M«=l-8, 
Re = 10 s , h = 0-3r (h and r are local strake height and cone radius, respectively.). 
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The side strakes push the primary vortices away from the leeward plane of geometric 
symmetry, and hence they prevent communication of the flow disturbance from the 
two sides. It should also be noticed that the C p curves of Figures 1 and 3 show that the 
side strakes provide additional lift besides their function of controlling the flow 
asymmetry. 

5.3. Unsteady Asymmetric Flow Around a Circular Cone 

Keeping the Mach number at 1-8 and Reynolds number at 10 s , the angle of attack is 
increased to 30° for the flow around the circular cone of Figure 1. Figure 4 shows the 
results of this case using the FDS scheme with the thin-layer NS equations, 
FVS-scheme with the thin-layer NS equations and the FDS scheme with the full NS 
equations. Here, we show the history of the residual error, the lift coefficient and the 
total-pressure-loss contours. For the first solution, pseudo-time stepping was used up to 



Thin-layer eqs. (FDS) 


x 1(T 3 



Thin-layer eqs. (FVS) 



Full Navier-Stokes eqs. (FDS) 

Figure 4. Comparison of unsteady asymmetric flow solutions with vortex shedding around a circular cone; 
periodic flow response; a = 30°, M* = T8, Re = 10\ At = 10 -3 . 
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10,000 iterations and the solution was monitored every 500 iterations. The solution 
showed that the asymmetry was changing from the left side to the right side, which 
indicated a possibility of unsteady asymmetric vortex shedding. The residual error was 
also oscillating. The computations were repeated starting from the 3,500 iteration step 
using time-accurate calculations with At = 10 -3 . The residual-error and lift-coefficient 
curves show the time history of the solution. It is seen that the residual error and the 
lift coefficient show a transient response which is followed by a periodic response. 
Snapshots of the total-pressure-loss contours are shown at time steps of n = 15,000; 
15,200; 15,400 and 15,700. At n = 15,000, the asymmetric flow is seen with an already 
shed vortex from the right side. As time progresses, the shed vortex is convected in the 
flow and the primary vortex on the left side stretches upwards, while the primary 
vortex on the right gets stronger and expands to the left side. At n = 15,700, the 
primary vortex from the left side is shed in the flow field. It should be noticed that the 
solution at n = 15,700 is exactly a mirror image of that at n = 15,000. The solution 
from 15,000-15,700 represents one half the cycle of shedding. The solution from 
15,700-16,400 represents the second one half of the cycle (not shown). The periodicity 
of the shedding motion is conclusively captured. The period of oscillation is 
10 3 x 1,400 steps = 1*4 which corresponds to a shedding frequency of 4*488. 

Figure 4 also shows the results of the FVS scheme with the thin-layer NS equations 
for one-half cycle of oscillation. Using the FVS scheme, the flux limiters are turned on 
initially and, as can be seen from the logarithmic-residual curve, the solution becomes 
symmetric and steady after 5,000 time steps. Next, the flux limiters are turned off, and 
the solution shows a transient response up to 12,000 time steps. Thereafter, the 
solution becomes periodic with periodic asymmetric vortex shedding. The solution was 
monitored every 100 time steps, and we show snapshots of the total-pressure-loss 
contours between n = 13,900 and n = 14,600. Although the process of adjusting the 
time instants in order to match those of the FDS solution is difficult, it is seen that the 
captured snapshots of the FVS solution almost match those of the FDS solution. 
Comparing the FVS solutions at n = 13,900 and 14,600, it is seen that they are mirror 
images of each other. Hence, periodic flow response has been achieved with a period 
of 1,400 x 10 3 = 1*4; which is exactly the same period of shedding as that of the FDS 
solution. This pinpoints the high numerical dissipation effect of the FVS scheme when 
the flux limiters are also turned on. The resulting numerical dissipation in the FVS 
scheme is large enough to dampen the random disturbances of the flow solution. By 
turning off the flux limiters in the FVS scheme, the random disturbances grow, 
producing the asymmetric unsteady vortex shedding. This also shows that the FDS 
scheme, even with the flux limiters turned on, is less dissipative than the FVS scheme. 
Finally, we show the results of the FDS-scheme with the full NS equations. The 
solution conclusively confirms the previous solutions and hence the unsteady solution is 
not scheme- or equation-dependent. 

Figure 5 shows snapshots of the surface-pressure coefficient using the FDS and FVS 
schemes with the thin-layer NS equations. They are in excellent agreement with each 
other. The oscillation of the side force is clearly seen in this figure. 


5.4. Unsteady Asymmetric Flow Around an Elliptic Section Cone, f r = 0*8 

Figure 6 shows the solution for an elliptic-section cone with fineness ratio of 0*8. The 
residual-error curve shows that the solution produces a symmetric flow through the first 
5,000 time steps. Afterwards, the solution shows a transient unsteady flow response for 
2,500 time steps which is followed by an unsteady, perfectly periodic, vortex-shedding 
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Figure 5 Comparison of snapshots of surface-pressure coefficients around a circular cone; a =30°, 

M„ = 1-8, Re = 10 s , A/ = 10 -3 . 


solution. The lift-coefficient curve shows the same nature of the solution as that of the 
residual-error curve. This case is carried out using time-accurate stepping with 
Af = 10 -3 . 

We also show snapshots of the total-pressure-loss contours and surface-pressure 
coefficients at the time steps of 12,000, 12,500, 13,000, 13,500, 14,000 and 14,500. The 
solutions at n = 12,000 and 14,500 are mirror images of each other which confirm that 
the solution is periodic. The period of oscillation is 5,000 x 10" 3 = 5 which corresponds 
to a shedding frequency of 1-257. At n = 12,000, the total-pressure-loss contours show 
that the right-side vortex is stretched, having two vortices; one at the top and the 
second one below it. In addition, a secondary vortex is seen at the surface. The 
left-side vortex has expanded to the right with two vortices beneath it. At n = 12,500, 
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Figure 6. Unsteady asymmetric flow solution with vortex shedding around an elliptic-section cone; 
periodic flow response; cr = 25°, M„— 1-5, Re — 10 5 , /, 0-8, At 10 


the top vortex on the right side has been almost shed while the one below it is 
expanding. At n = 13,000, the top vortex on the right side has been shed and 
converted with the flow, while the vortex below it is expanding to the left. As time 
passes, the vortex on the left side is stretching upwards and the vortex on the right side 
is expanding to the left, as seen from the snapshots at 13,000, 13,500 and 14,000. At 
n = 14 500 the vortices on the left side and right side become mirror images of those 
on the’ right side and left side, respectively, at n = 12,000. Again, this case conclusively 
shows that unsteady vortex-shedding flows are captured. 
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5.5. Unsteady Asymmetric Flow Around an Elliptic-Section Cone, f r 0-2 
(Thick Wing Case) 

This case is presented to show that asymmetric vortex shedding also > exists Tor ^ng-like 

sections. Here, the elliptic-section fineness ratio is reduced to 0-2. To obtain this 

impressive flow case, we have to decrease the free-stream Mach number to 14, 
imprest! vc uu , , ;nrr _ ase the free-stream Reynolds number to 

2 ^r^r&justmems have been made to enhance the critical Bow conditions for 

flow asymmetry. For the same section flneness ratio, same angle of attack and same 

riel stream Mach number and for free-stream Reynolds number range of 10 5 - 
tree stream m etcadv asvmmetric. It is unsteady only when 

1-5 x 10 6 , the flow is either symmetric or steady asymmeun-. , j 

the Reynolds number is increased to 2 x 10 6 . 
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Figure 7 shows the time-accurate (A t = 2 x 10 3 ) results of this case w ic inc u es 
the time-history of residual error, the time-history of the lift coe cient, snap s o so 

the total-pressure-loss contours and snapshots of the SUI ^ e *P^ es ^ r ^i oe Cien *' 6 

snapshots at n = 15,000, 15,100, 15,200, 15,300, 15,400 and 15,500 repres p- 
proximately one-half the cycle of the periodic flow response. e o a 
contour at n = 15,000 shows that the left-side vortex « stretched while the nght-side 
vortex has expanded, covering a large region of the left si e o ® ov ^ » 

the wing. Under the right-side vortex, a strong secon ary vo ex is • 

n = 15, 100, the left-side vortex shows two regions of vortical flows; one a p 

another one below it. Both vortex regions of the left vortex rotate m the same 
clockwise direction. At n = 15,200, the top vortex is shed into the fi eM, *™ s h e d 
one below it gets stronger and stretches upwards. At rt , ’ ’ h 

vortex from the left side is converted in the #ot», the left vortex '^xpand.ng to the 
right and convecting vorticity to the nght vortex. The right vortex l g g g ’ 
shrinking in thickness and stretching upwards. A secondary vortex is forming u 
left vortex, and the secondary vortex under the right vortex « *mm.sh nng, a 
n = 15, 100, the flow is approximately a mirror image of that at n - 15,500. 1 he num 
of time steps for one cycle of periodic response is 1050, which giwss a period o 
oscillation of 2 x 10 -3 x 1,050 = 2-1 corresponding to a frequency of 2-992. 

5.6. Unsteady Asymmetric Flow Around a Diamond-Section Cone, / 4 = 0-2 (Thick 
Wing Case) 

Figure 8 shows the results of the time-accurate (Ar = 5 x 1(T 4 ) results for this case 
which include snapshots of the total-pressure-loss contours and snapshots of the 
surface-pressure coefficients. The snapshots are given at n = 11,500, 12,000, 12,500, 
13 000 13 500 and 14,000. The number of time steps for one cycle of periodic response 
!s5, which gives a period of oscillation of 5 x 1<T 4 x 4,500 = 2-25 correspond, ngto 
a frequency of 2-793. It should be noted here that the angle of attack is 38 , which is 
higher than that of the elliptic-section cone of Figure 7, where the angle of attack is 34 
and all the other flow conditions are the same. Comparing the surface pressure curves 
of the elliptic-section wing (Figure 7) and the diamond-section wing (Figure), we 
conclude that the diamond-section wing has less asymmetric strength and higher lift 
coefficient than those of the elliptic-section wing. 

6. CONCLUDING REMARKS 

The unsteady, compressible, thin-layer, Navier-Stokes equations are used along with 
several computational schemes to numerically simulate steady and unsteady asym- 
metric vortex flows around cones. The concept of the locally conical flow assumption 
has been developed and discussed. A steady asymmetric flow solutions has been 
presented and verified for a circular cone. Passive control of the flow asymmetry has 
been demonstrated for the circular cone by using side strakes. Unsteady, asymmetric 
vortex flows with periodic vortex shedding have been presented for cones with a 
circular section, an elliptic section of fineness ratio of 0-8, an elliptic section of fineness 
ratio of 0-2 and a diamond section of fineness ratio of 0-2. The unsteady asymmetric 
flow solution for the circular cone has been verified using two schemes with the 
thin-layer and full NS equations. The present study shows that, for the same flow 
conditions, circular section cones produce the strongest flow asymmetry while the 
diamond section cones produce the weakest flow asymmetry. It is conclusive that 
unsteady flow asymmetry with vortex shedding has been captured. It should be noted 
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that these solutions are based on the locally conical flow assumption and hence they 
must not be used for quantitative comparisons with three-dimensional flow results. The 
reason behind such a restriction is the length scale involved with the unsteady and 
viscous terms of the locally conical NS equations. However, the solutions are 
computationally economical for qualitative and parametric studies of asymmetric flows. 
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The unsteady, compressible, thin-layer Navier-Stokes equations are used to numencally study 
the passive control of steady and unsteady supersonic asymmetric flows around circular and 
noncircular cones. The main computational scheme of the present study is an implicit upwind, 
flux-difference splitting, finite-volume scheme. Passive control of flow asymmetry is studied by 
using a vertical fin in the leeward plane of geometric symmetry and side strakes with and wi ou 
thickness at different orientations. The study focuses on circular-section cones since they are the 
most likely section-shapes for strong flow asymmetry. Side-strake passive control is shown 
more efficient and practical than vertical-fin passive control. • .992 a^c Pm*. Inc. 


Introduction 

In the moderate to high angle-of-attack (AOA) ranges, which are typical 
flight conditions for highly maneuverable fighter aircraft and missiles, extensive 
regions of vortex-dominated flow develop on the vehicle. 

Within these AOA ranges, the cross-flow velocity components and the gra- 
dients of other flow variables become of the same order of magnitude as or 
higher than those of the axial direction. Consequently, flow separation occurs 
and vortices emanate from the three-dimensional separation lines of boundary- 
layer flows on wings, strakes, and fuselage of the vehicle. If the vortices i are 
symmetric and stable, their influence can be exploited favorably to pro 
high lift and maneuverability for the vehicle. On the other hand, if the vortices 
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become asymmetric or if vortex ^mttricTftog' fora. and 

corresponding ya“oiung!and pitching “^*1 

— n is another un- 

favorable vortexondnced phenomenorr ^ ^ and pointed slender 

Highly swept, round- and sharp- of fichter aircraft and missiles, 

bodies are common aerodynamic these isolate d aerodynamic com- 

The study of vortex-dominated flow vortex-dominated flows. The so- 

ponents adds to our basic understan 1 8 y^cs the high AOA rang e 

lution of asymmetric vort * x OW a controllability of fighter aircraft and mis- 
is vital to the dynamic stabi 1 y a h relative incidence (ratio 

siles. The onset of flow exceeds certain 

of angle of attack to nose sem pe & relative incidence, flow asymmetry 
critical The origin of natural 

develops due to nalu ” j ;iWli a n acoustic disturbance, or a similar 
disturbances may be a trans disturbances, however, are caused by 

disturbance of short durat . • the nose geometric symmetry 

tion shape are ™P 0 ^ M /“^ve^nmade to simulate asymmetric 
Several computabonal attempts have been m ^ ^ (J] ^ 

vortical flows around slender bodies [ J’ along ^th two different im- 

the unsteady thin-layer Navier- o around cones with different 

phci. schemes to simulate t = * ^“as “Ised on a 5«-semi- 

cross-section shapes. Th ^ conica i flow assumption. The first 

apex angle circular cone difference splitting, finite-volume 

computational scheme was an finite-volume 

scheme and the second one wa number being held constant at 1 .8 

scheme The Mach number .0° to 30°. At 

and 10 5 ; wL obtained and the results of the two 

a - 10 , a steady sy * . _ 70 ° and irrespective of the type 

schemes were in excellent agteem - so | ul | on was obtained and 

or level of the disturbance, a steady asy ™"V Two types of flow 

the results of the two schemes were ““^"or a random tranca, ion- 
disturbances were used; a random r sideslip disturbance with short 

error disturbance, and a ^rolled —t si 'SmrbtTthe solution was 
duration. For the controlled rans. ^ent sideriip < Jisturh was als0 

unique, and for the uncontrolled — sides on 

unique with * h i" n un t e' d iymmetric solution with vortex shedding 
was^btained and the vortex sh Jding was perfectly periodic. Next, the angle 
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. pnt fixed at 20° and the Mach number was increased from 1 . 
of attack was k P The solutions showed that the asymmetry became 

C"" a.so presented for cones with elliptic and diamond 

cross-sectional areas. 161 several issues related to the 

a later paper by the P res \"‘ f^/.ifshown that a unique asym- 

asymmetric flow solutions were addressed^ minimum grid 

metric flow solution is obmmediW^of ~ ^ ^ ^ the 
spacing at the solid boundary. ^ shown that for the same flow 

random nature of the disturbance. It was ^ wj(h sharp 

conditions and same section finen - CQnes Agam . it was also 

cones were presented and audmd. been mad e to control asymmetric 

Experimental research efforts h CVTTimetr ic forces and the resulting 

flows for eliminating or attenuatmg the |vMontro , [11-13] 

moments by using either P— ° of a vertic al fin on the 
methods. Passive-control me 0 symmetry [8], the use of fixed or 

leeward side along the ^plane o g D f a rotatable forebody tip having 

movable forebody strakes [9 2], its base to an elliptic shape at 

variable cross section (from a circuia ^ inc l ud e the use of blowing ports 
its tip [10]). Active control th^forebody surface [11, 12]. 

both passive [5] and active contro pulat ional study for passive control 

In the present paper, we present ^ ^ ^ usjng a ver tical fin in 

of supersonic asymmet tr v and side strakes with and without 

the leeward plane of geometnc symmeUy and side^ ^ ^ ^ focused 

^"mular-secuon cones, a few applications have been 

^nciripred for noncircular-section cones. 


Formulation 


Governing Equations . 

The three-dimensional, 

^"ations’in"' Ind” dfre^ons yields the transformed equations 
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do s£, X&h _ djgjs = n- 1,2,3, (1) 

Tt ae 


where 

>* 1 

Q = - = - [p, puu P u 2, pu 3 , pe]' ( 2 ) 

J J 

E s = inviscid flux 

= - [pl/„ PUiUs + dtfp, pu 2 U s + d 2 ?p, pu 2 U s + d&p, (pe + p)U s f, 

s = 1, 2, 3 (3) 

(£ v ) 2 = viscous and heat-conduction flux in the £ 2 direction 


— 

= - [0, d k £r k 1 , d k £ 2 T k2 , d k £ 2 T ki , d k £ 2 (u n T kn Q k )] 

J 

(E v ) 3 ** viscous and heat-conduction flux in the | 3 direction 

(4) 

— 

= - [0, d k ^r k u d k t?T k2 , dic^Tw, d k £ 3 (u„ T kn ~ Qk)] 

(5) 

— 

U s = d k eu k . 

( 6 ) 


The first of the three momentum elements of Eq. (5) is given by 



1 A/qcM / 1 n v3 t i 

dti Tk ' " 'rT V 't *W)' 

(7) 

- 

where 

1 a du k 

<t> = dk?d&\ ^=3 d £ ^3- 

( 8 ) 


The second and third 
subscript 1 , everywhere 
of Eq. (5) is given by 


momentum elements are obtained by replacing the 
in Eq. (7), with 2 and 3, respectively. The last element 


* 


d k t\u„Tkn ~ Qk ) 

A/f r n s ,7 1 

- ^ + *[2 ap (M? + ^ + ^ + (7 - DPr ^ 3 . 


(9) 


For Eq. (4), in the case 
the elements are given 


of thin-layer approximations applied to the £ 2 direction, 
by equations similar to Eqs. (7H») th< = exce P tl0n 
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• ,3 u >-2 T u; c c Pt of thin-layer Navier-Stokes equations is used 

of replacing £ by £ . cases using a ver tical fin since the existence of 

only for the passive con 1 ^ ^ which perpe ndicular to the cone’s 

the fin creates a second ^ g directlon . For the passive 

thin layer. This viscous ay equations in £ 2 direction are used 

control cases using side slrakes ' ' J j the 52 auction due to the method 
since the viscous layers on the strokes are mt s parameters for 

used to construct the gnd around the brakes. 1 ne ^ ^ for , he 

the dimensionless form oftheequa ion respectively. The 

length, velocity, ume density, and the pressure, p, is related 

( 10 ) 


p = {y- l)p^-^(w? + «2 + M 2 3)j. 

The viscosity is calculated from the Sutherland law 

>/ l t£ \ C = 0.4317, 
\T+C/’ 


_ J"3/2J 


(ID 


rrssK ssf* 

-^jr,s=s,t5aM53 « - «»» 


Boundary and Initial Conditions 

Boundary conditions are explicitly '^“nTlt Sane" of geome^c 

outflow conditions and solid- oun ary At the f ar .fleld inflow boundaries, 
symmetry, periodic conditions are used fM the faHd ^ ^ rf 

f r St ?r n^AUheX S Sfl ot boundaries first-order extrapolation from 
the solution. At the tar boundary, the no-slip and no-pene- 

the intenor points is used. On = = 0) and the normal pressure 

tration conditions are enforce ( i 2 ature t h e adiabatic boundary 

to uniform flow with u^ u 2 3 boundary conditions are enforced 

For the passive control applications, solid-boundary conui 

on both sides of the fin or the strakes. 

HIGHLIGHTS OF THE COMPUTATIONAL SCHEME 

anTmp^h "upwi^^ux^ffe^M spfi^ng^ ^ite-volume s^^^f crnpl®^ 
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the flux-difference splitting scheme of Roe, which is based on the solution of 
the approximate Riemann problem. The inviscid flux difference at the cell 
interface is split into left and right flux differences according to the signs of 
the eigenvalues of the Roe averaged-Jacobian matrices of the inviscid fluxes 
A = dE/dq ■ 5=1-3 Flux limiters are used to eliminate oscillations in the 
shock region. The viscous and heat-flux terms are centrally differenced. The 
resulting difference equation is solved using approximate factorization m the 

£ 1 2 ? an( j p3 directions. The computational scheme is coded in the computer 
program “CFL3D ” In this program, an implicit, flux-vector splitting, finite- 
volume scheme, which is based on the Van Leer scheme [14], can also be 
used instead of the flux-difference splitting scheme. The flux-vector splitting 
scheme is also used to solve for the unsteady asymmetric flow application in 
this paper. This application is a validation of the solution obtained previously 
[5] for the same application using the flux-difference splitting scheme. 

Since the applications in this paper cover conical flows only, the three- 
dimensional scheme is used to solve for locally conical flows. This is achieved 
by forcing the conserved components of the flow vector field to be equal at 
the two axial planes located at =0.95 and 1.0 of the conical grid. 


Computational Studies 

1. Passive Control For a Circular Cone Using a Vertical Fin, a = 20° 

In this section, we consider the control of steady asymmetric flow around 
a 5° -semiapex angle circular cone at an angle of attack a — 20 , freestream 
Mach number M . = 1.8, and freestream Reynolds number Re = 10 s . Two 
vertical fins of heights h = 0.5r and r are placed in the leeward plane of 
geometric symmetry, where r is the cone local radius. Figure 1 shows the 
results of this study, which include the total-pressure-loss contours, the surface- 
pressure coefficient versus the angle 6 ( 6 is measured from the leeward plane 
of geometric symmetry in the clockwise direction), and the lift coefficient 
versus the number of iterations. The figure also shows the results of the asym- 
metric flow without a vertical fin [5]. With h = 0.5 r, the flow is still strongly 
asymmetric after 34,000 iterations. Comparing the case of no fin with the 
case of h = 0.5r, it is seen that the asymmetry changed sides. This is due to 
the random nature of the disturbance, which is a machine round-off error. 
With h = 0.5r, two vortex cores, which are connected to each other and to 
the body through free-shear layers, develop from the left side of the body. 
From the right side of the body, a free-shear layer develops and crosses over 
the fin to the left side of the body. It produces two vortex cores; one at each 
comer of the body-fin juncture with secondary separations below them. This 
case has been solved accurately in time but it does not show any vortex shed- 
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Iter otions 




iterations 




tier otions 


Fig 1 Passive control of asymmetric flow around a circular cone using a vertical fin. a 20 

M = 18? Re = 10 s , h = 0.5r and r (r = cone local radius). 


dine or unsteadiness. When the fin height is increased to h - r, perfect flow 
symmetry is obtained. The lift-coefficient curves show that when flow asym- 
metry develops, the lift coefficient increases over a small number of iterations 
an^remains constant thereafter. When the flow becomes symmetric, as with 
h = r the lift coefficient does not increase. The reason ^hindthe flow 
me try with h = 0.5 r is that the free-shear layer from the right-hand side of 
th^body is still higher than the fin height, which allows the flow disturbances 

from the right and left side to interact. 

Fieure 2 shows a typical grid for studying passive control using a vertical 

fin It contains 161 X 81 grid points in the wrap-around and normal directions, 

respect^ The grid is generated by using a modified Joukowski transfor- 
mation with clustering in the normal direction at the cone surface, an c u 
tering hvthe wrap-around direction at the venicai fin surfaces. The maxtmutn 
Sins of the computational domain is 21r. The figure allows " 
of the cross-flow velocity in the cone-fin-juncture region for h - 0.5 r and r. 
It is noticed that two small recirculating bubbles exist under the vortex cores. 
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, Passive Control For a Circular Cone Vsing Side Strokes, a = 20« 

In Fig. 3, control of flow h = 

irst application is consl ^ e flo w perfectly symmetric since the two primary 
^ preventing the flow dis- 

vortex cores, lett ana ngni, <u v seen t hat there are four 

turbances of the two sides from m er * three are secondary vortices, 

vortices on each side; one is a primary v ° 0 d 2 70°, where the 

The pressure-coefficient curve *ows jumps a. 6 90 MO , , ^ ^ 

strakes are ! ! ^Sn-SnDo. case. Moreover, the lift coefficient 

^r n s,d«. ln.r ol case is double that of the vertical-fin-control case. 
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1*i io-» 
12 - 




FIG. 3. Passive control of asymmetric flow around a circular cone using sharp-edged thick 
strakes. a = 20°, M w = 1.8, Re = 10 5 , h = 0.3r. 


It is concluded that side-strake control in comparison with vertical-fin contro 
not only is efficient in providing higher lift, but also is more practical due to 
the strake’s shorter height, and hence lesser weight. 

3. Passive Control For a Circular Cone Using Side Strakes with and 

without Thickness, a = 30° 

In this case the cone angle of attack is increased to 30° keeping all the 
other flow conditions fixed. This flow application has been solxed previously 
in Ref [5] by Kandil, Wong, and Uu using a flux-difference splitting (FD ) 
scheme The results showed unsteady asymmetnc flow with periodic vorte 
Ihedchng The total-pressure-loss contours of the time steps from 15,000 to 
^ 700 ™ representing one-half the cycle during the periodic response, is shown 
!„ ft 4 <^er unsteady asymmetric flows with periodic vortex shedding 

around elliptic-section and diamond-section cones have a '”^ pr “ was 
hv the authors in Refs. [5, 7]. In all these applications, the FDS scheme 
used In order to show that the unsteady asymmetric solutions are not scheme 
dependent^the effect of computational methodologies and numenod 
pation on the solutions are examined in the present paperThec^ of asym 
metric flow around a circular cone has also been computed using the flux 








PASSIVE CONTROL OF FLOWS AROUND CONES 


89 


1 i y~~ 


o *« 

—t 

-8 

-«0 



— to — is ao* 5110 ' 

Iterotions 
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m — Is k n , to* 



15,500 


15,600 


15,700 


Iterations 15,400 

„ t - r flow with vortex shedding around a circular cone using FDS 

rrr -o- 

vector splitting (FVS) scheme rflhe same ^ £££ m 

mic-restdual curve, the so u l tumed off and the solution shows a 

time steps. Next, the flux li Thereafter, the solution becomes 

transient response up to 12, shedding. The solution is shown 

periodic with periodic asymm ^ un til time step 14,600. 

every 100 time steps srarung from time riep^3,0W u ^ ^ ^ ^ 

Although the process ofedjustingi ‘ 0 f the FVS solution 

solution is difficult, st^ps .5,000, .5,100, .5,200 

almost match those of the <- 7ft0 res nectively. Comparing the FVS 

15 ,300, .5,400, it is seen that they ate 

solutions at time step of 1 3.9W flow responS e ha s been achieved 

mirror images of each oth .H^ = . P wh , ch is exactly the same period of 

with a period of 1 ^00 X it ) ^ pinpoints the high numencal dis- 

d :rno n w m soM,on. By turning off the flux limiters in 
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F,c. 5. Unsteady asymmetric flow with vottes shedding aronnd a circular cone using FVS 
schemes, a = 30°, M „ = L8, Re - 10 s , A t 10 . 

the FVS scheme the random disturbances can grow producing the asymmetric 
S,-; shedding. This also shows , ha, the ^"^enw^h 
the flux limiters turned on, is less dissipative than the FVS scheme, inese 
results conclusively show that unsteady asymmetric flows are obtained irre- 

^Nex" we^nsidCT^^cmttd^^^unsteady asymmetric flow using sha*- 

edged thick strakes and of this 

the strake shapes, the height is kept at 0.3r. Figure e > sn o 
,tudv using sharp-edged and flat-plate strakes at 6 = 0 ,10 ,and 1U , wnere 
5 is the angle measured in the counter-clockwise direction from the honzonta 
ine at 0 = 90°. For this angle of attack, all the strake orientations are still 
effective in eliminating the unsteady asymmetric vortex shedding and ren- 

dering the flow P e ^ 1 Y all controlled 

b 7 SMtanTha, Of the asymmetric flow case. Wi.h si.gh, 
differences in the lift coefficient, the highest lift is produced by the flat-plaK 
strakes with 5 = -10°, where the primary vortex cores are sligh y c 
the body surface than for the other cases. 
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Comparing the results of the sharpredged thick strakes at a = 30 with 
those of the sharp-edged thick strakes at « = 20° (Fig 3), we notice that the 
primary vortex cores of the former are closer to the plane of symmetry and 

hieher above the cone surface than those of the latter. 

fn the bottom row of Fig. 7, we show typical gnds for the cases ofsharp- 
edeed thick strakes and the flat-plate strakes with 8 - 10 . The gnds are 
generated by using a hyperbolic grid generator with transfinite gnd mterpo- 
lation to refine the grid in the strake region. 

4. Passive Control For a Circular Cone Using Sharp-Edged Thick Strakes, 
a = 40° 

For this case, the angle of attack is increased to a = 40° keeping all the 
other flow conditions fixed. The same sharp-edged thick strakes ofthepre ™^ 
case (a = 30°) have been used along with the same gnd. Figure 8 shows the 
Sts of this case. It is seen that although the C p -curve looks perfectly sym- 
metric and although the lift coefficient curve does not show any increase after 
the time step 4,000, the total-pressure-loss contours show very^ slight ^ym- 
metrv. This indicates that the current height of the strakes migh no 
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ficient to yield How symmetry at higher angles of attack. The present solution 
*ows that the vortical flow substantially stretches upwards. 

5 . passive Control For an Elliptic-Section Cone Ustng a Vertical Fin, 

= 25° 

, . nrof i U ce a substantial flow asymmetry (of the same order as 

In order to p . r an ellintic-section cone of fineness ratio 

that of the circular cone of Fig. 1) P _ anc i t ^ e f re estream 

fr = 0.6, the angle of attack has been mcreasedto «-25 flow 

Mach number has been decreed to , M - 1.5. ■ c 2a is the 

“YT 

rotatd no^ecbons, and was generat^ by a m^ loukowsk. 

transformation with m ' n “ U ” *owi in Figs. 9 and 10. For the 

(for the shows an mcrease near step 7,000, and it 
case with no fin, t = { 5a the flow is still strongly 

ramains constant increase near step 8000. The 

asymmetric and the show a very i ong free-shear layer from the left 

total-pressure-loss c free . shear layer becomes higher than the vertical 

f "rr^otr the fin to the left side. Two primary-vortex cores are 



No Fin 


h=1.5a 


h=2a 






F,G. 9. Passive control ofasymmetnc minor-axis length). 

= 25°, Moo = 1-5, Re = 10 , tr - v.o, 
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F,G. I OXompamon of I, ft ^^‘^7 


satssr- 

sides, and flow symmetry will be achieved. 

6. fturiw Comro/ For a Diamond-Section Cone Vsmg a Vert, cal . 

For this case, the section i^Ts' With the exception of the section 

25°, and the fteestream Mach number K ^ ^ ^ ^ of ^ 

fineness ratio, the flow conditions section fineness 

elliptic-section cone. This simply shows JWfcrt nes with shatp edges pro- 
mt. o and same flow condmon^ dmmo ^ion cones, 
duce less flow-asymmetry strength ^ cone flow with and 

Figure 1 1 -f” e r S “he fljw^omrol case, a symmetric flow has been 

without a vertical fin. , • ht l = \ ^ a , which is shorter than that 

obtained using a vertical fin of height h za, 

«™;rprt for the elliptic-section case. 


Concluding Remarks 

onntrni of steadv and unsteady, super- 
Computational studies for passive Qu , using vert i ca i fins in 

£ts£ ssr ■ 

rr “ ^ *- *- — uang 
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Fio 11 Passive control of asymmetric flow around a diamond-section cone 
fine's" 25°, = 1 .5, Re = 10 s , fr = 0.8, h = 


using a vertical 


an implicit, upwind, 

control Symmetry. It has been shown that 

* — " ,h the 

irarilCdl fin 
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